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ABSTRACT 

We present global multi-dimensional numerical simulations of the plasma that per- 
vades the dark matter halos of clusters, groups, and massive galaxies (the 'intracluster 
medium;' ICM). Observations of clusters and groups imply that such halos are roughly 
in global thermal equilibrium, with heating balancing cooling when averaged over suffi- 
ciently long time- and length-scales; the ICM is, however, very likely to be locally ther- 
mally unstable. Using simple observationally-motivated heating prescriptions, we show 
that local thermal instability (TI) can produce a multi-phase medium — with ~ 10^ K 
cold filaments condensing out of the hot ICM — only when the ratio of the TI timescale 
in the hot plasma (txi) to the free-fall timescale (tg) satisfies ^Ti/iff ^ 10. This criterion 
quantitatively explains why cold gas and star formation are preferentially observed in 
low-entropy clusters and groups. In addition, the interplay among heating, cooling, and 
TI reduces the net cooling rate and the mass accretion rate at small radii by factors of 
~ 100 relative to cooling-flow models. This dramatic reduction is in line with observa- 
tions. The feedback efficiency required to prevent a cooling-flow is ^ 10"^ for clusters 
and decreases for lower mass halos; supernova heating may be energetically sufficient 
to balance cooling in galactic halos. We further argue that the ICM self-adjusts so that 
tTi/tff ^ 10 at all radii. When this criterion is not satisfied, cold filaments condense out 
of the hot phase and reduce the density of the ICM. These cold filaments can power 
the black hole and/or stellar feedback required for global thermal balance, which drives 
tTi/tg ^ 10. In comparison to clusters, groups have central cores with lower densities 
and larger radii. This can account for the deviations from self-similarity in the X-ray 
luminosity-temperature {Lx — Tx) relation. The high- velocity clouds observed in the 
Galactic halo can be due to local TI producing multi-phase gas close to the virial radius 
if the density of the hot plasma in the Galactic halo is ^ 10~^ cm~^ at large radii. 

Key words: galaxies: clusters: intracluster medium; galaxies: halos. 



1 INTRODUCTION 

The cooling of baryons in dark matter halos is a critical driver 
of galaxy formation on large scales (e.g., Hoylc 1953; Rccs & 
Ostrikcr 1977; White & Rccs 1978). If the cooling time at the 
virial radius is shorter than the free-fall time, most of the gas 
does not virialize at large radii and instead remains cold and 
flows into the halo along cosmological filaments (Birnboim & 
Dckcl 2003; Kcrcs ct al. 2005). For halo masses larger than 
~ 10^^ M0 (roughly the halo mass of the Milky Way; Xue et 
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al. 2008), however, most of the baryons undergo a virial shock 
and produce a pressure-supported intracluster medium (ICM; 
for simplicity, we refer to the viriahzed plasma at all halo 
masses as the ICM). 

Even when cooling is negligible near the virial radius, it 
can become important at smaller radii where the gas den- 
sity is higher. Indeed, the cooling time of the hot plasma at 
the centers of groups and clusters is often observed to be 
much shorter than the age of the system (e.g., Cavagnolo ot 
al. 2008). In the absence of heating, the ICM would lose 
pressure support on a cooling timescale, leading to a mas- 
sive cooling-flow (e.g., Fabian 1994). However, the absence 
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of lower energy X-ray lines (e.g., Tamura et al. 2001; Molcndi 
& Pizzolato 2001; Peterson ct al. 2003) and the lack of prodi- 
gious star formation in cluster cores (e.g., O'Dea et al. 2008; 
Raffcrty, McNaniara, & Nulscn 2008) suggest that cooling 
is approximately balanced by some form of heating, at least 
when averaged over a cooling timescale. Heating is also re- 
quired to explain the high mass cutoff in the galaxy mass 
function relative to the dark matter halo mass function (e.g., 
Croton ct al. 2006; Bower ct al. 2006). 

Energetically, there are a number of heating sources that 
can in principle balance cooling in hot halos: e.g., thermal 
conduction in clusters (e.g., Zakamska & Narayan 2003), the 
kinetic energy of galaxies and mini-halos (e.g., Kim, El-Zant, 
& Kamionkowski 2005; Dckcl & Birnboim 2008), and the 
energy produced by a central active galactic nucleus (AGN) 
in the form of jets and bubbles of relativistic plasma (e.g., 
Birzan ct al. 2004). However, most non- feedback mechanisms 
are globally thermally unstable (e.g., Kim & Narayan 2003; 
Parrish, Quataert, & Sharma 2009). That is, even if heating 
balances cooling at a given time, small perturbations about 
that equilibrium will cause run-away heating/cooling of the 
core. AGN feedback, however, can be globally thermallj stable 
if the efficiency of feedback is sufficiently large (e.g., Guo ct 
al. 2008). Physically, this is because the accretion rate onto 
the central black hole increases as the cooling rate in the 
halo increases, providing the necessary feedback loop. It is 
plausible that supernova heating can provide a similar source 
of globally stable heating in galactic halos with ongoing star 
formation. 

Although AGN feedback can plausibly ensure global ther- 
mal stability of cluster cores, the ICM is very likely to be lo- 
cally thermally unstable because of the strong density depen- 
dence of the cooling rate (Field 1965). The situation is analo- 
gous to the hot phase of the interstellar medium (ISM) which 
is thermally unstable but is maintained by supernova heating 
(McKec & Ostrikcr 1977). While the thermal instability (TI) 
has been invoked to explain a number of observations of cool- 
core clusters in the past, previous treatments have typically 
focused on models without heating (e.g., Nulscn 1986; Kauf- 
mann ct al. 2009). Other studies, very similar in spirit to 
ours, have emphasized the importance of AGN feedback due 
to accretion of cold gas, but have focused on the survival 
of dense, cool blobs once they form, rather than their origin 
(e.g., Pizzolato & Soker 2005). Here we focus on how local 
TI produces cold filaments condensing out of the hot phase 
for virialized halos in rough thermal balance. 

Observationally there is a sharp threshold in central en- 
tropy/cooling time below which cold gas is observed in groups 
and clusters (e.g., Cavagnolo ct al. 2008; Hcckman ct al. 
1989; Hu, Cowie, & Wang 1985; McDonald ct al. 2010). 
These cool-core clusters also preferentially show evidence for 
AGN feedback and star formation (e.g., Rafferty, McNamara, 
& Nulscn 2008; Hicks, Mushotzky, & Donahue 2010). The 
strong correlation between cold gas indicators and the prop- 
erties of the hot ICM motivates a model in which TI is the 
source of much of the cold gas in groups and clusters. The con- 
densation of cold gas in hot halos has also been invoked (e.g.. 



Mailer & Bullock 2004; Sommcr-Larscn 2006; Peck, Put- 
man, & Sommcr-Larscn 2008) as a source of the high velocity 
clouds (HVCs) of neutral hydrogen observed in the Galactic 
halo (e.g., Wakker & van Wocrdcn 1997) and in the halo of 
M31 (Thilkcr et al. 2004). 

In a previous paper (McCourt ct al. 2011; hereafter Pa- 
per I) we studied the physics of local TI in gravitationally 
stratified plasmas, including observationally-motivated heat- 
ing prescriptions, anisotropic thermal conduction, and mag- 
netic fields. We used local Cartesian numerical simulations to 
isolate the key physics. The unusual feature of the TI in the 
ICM (relative to, say, the molecular/atomic phase of the ISM) 
is that the growth time of the TI in the hot plasma {tTi) can 
be comparable to, or longer than, the free-fall time in the clus- 
ter potential (iff). In Paper I, we showed that the ratio of these 
two key timescales determines whether cold gas can condense 
out of the hot ICM. If the cooling time is short (txi < iff), ini- 
tially small density perturbations become nonlinear and cool 
blobs reach the thermally stable state at ~ 10^ K. These dense 
gas-rich filaments then fall through the hot ICM at roughly 
the free-fall rate (as in Pizzolato & Soker 2005); they may 
cool further, forming stars, and/or powering a central AGN. 
On the other hand, if the TI growth time is longer than the 
free-fall time {Iti > iff), the TI saturates before the pertur- 
bations become nonlinear and the TI cannot produce a multi- 
phase medium starting from small initial perturbations. Note 
that this criterion for TI to produce a multi-phase medium is 
analogous to the criterion for cosmological accretion to pro- 
ceed via cold streams rather than virialized plasma. However, 
our TI criterion applies at all radii in hot halos, not just at 
the virial shock. 

In this paper, we extend the calculations of Paper I to 
more realistic (but still non-cosmological) global models of 
the ICM in dark matter halos. The most uncertain input into 
these calculations is our treatment of the heating of the ICM. 
Observations of groups and clusters imply that there is an 
approximate balance between heating and cooling in clus- 
ter cores. Motivated by this, we explore two different heat- 
ing models; both assume that the heating is isotropic, thus 
sidestepping the very important but difficult problem of how 
energy is redistributed throughout the ICM. Our first model 
determines the heating rate from the average cooling rate at 
each radius in the halo (see also Paper I). We refer to these 
simulations as our 'idealized heating' simulations. Our second 
model assumes that the heating of the ICM is proportional to 
the mass accretion rate at small radii. We refer to these sim- 
ulations as our 'feedback' heating simulations (though they 
are also very idealized!). Fortunately, these two models give 
similar results. This implies that there are robust properties 
of hot halos that arise primarily from the interplay between 
local TI and global thermal equilibrium, and which are not 
very sensitive to the details of how global thermal equilibrium 
is maintained. 

Our calculations do not include the physics necessary to 
model the detailed properties of the gas-rich filaments pro- 
duced by TI. In reality magnetic fields, cosmic rays, and ther- 
mal conduction along magnetic field lines strongly infiuence 
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the geometry and energy balance in filaments (e.g., Fcrland ct 
al. 2009; Sharma, Parrish, & Quatacrt 2010). Star formation 
needs to be incorporated as an excitation mechanism for the 
observed Hq emission in at least some filaments (e.g., Mc- 
Donald ct al. 2011). Moreover, the physics of turbulent heat- 
ing/mixing and non-radiative cooling are probably required 
to explain the absence of some of the lower energy X-ray 
lines (Werner et al. 2010; Sanders et al. 2010); this physics is 
not well-resolved in our simulations. A detailed study of these 
processes is important, but beyond the scope of the present 
paper. Nonetheless, the simulations presented in this paper 
do accurately describe the conditions under which gas-rich 
filaments form out of the hot ICM; they also provide insight 
into the role of TI in determining the observed density and 
entropy profiles of hot halos, and in maintaining such halos 
in global thermal equilibrium. 

In Paper I we showed that magnetic fields and anisotropic 
thermal conduction do not significantly alter the evolution of 
the TI for halos in thermal balance. More specifically, conduc- 
tion changes the morphology of the resulting cold gas (turning 
blobs into filaments along magnetic field lines), but not the 
amount of cold gas produced by TI. As a result, we do not 
include magnetic fields and conduction in this paper. 

The interplay between local TI and gravity in a globally 
stable cluster cannot be captured in a one-dimensional model 
because the physics of cool overdense blobs falling in a hot 
pressure-supported atmosphere is absent in 1-D. Moreover, 
maintaining the pressure support of the hot ICM is critical to 
preventing a massive cooling-flow. The TI also evolves very 
differently in the presence of heating than it does in cooling- 
flows (Balbus & Soker 1989; see also Appendix A where we 
discuss the formation of multi-phase structure in cooling-flows 
in detail). These considerations imply that numerical stud- 
ies of hot halos should use multidimensional simulations and 
should include heating. Such simulations are the focus of this 
paper. 

The remainder of this paper is organized as follows. We 
describe the properties of our multi-dimensional simulations 
in section 2. Section 3 presents the results of our 'idealized' 
heating simulations. Simulations with simplified models of 
'feedback' heating are presented in section 4. We discuss im- 
portant astrophysical implications of our work in section 5. 
We conclude with a detailed summary of our paper in section 
6. Readers not interested in simulation details may skip to 
section 5. 



2 SIMULATION SETUP 

In this section we detail the equations that we solve numer- 
ically, and the initial and boundary conditions. We also de- 
scribe the cooling and heating prescriptions that we use. 



2.1 Governing Equations 

We solve the standard hydrodynamic Euler equations with 
heating, cooling, and external gravity: 



dp 
IE 



= -pV • V, 



dv „ 

p-^ = -^p + pg, 



e-ln (p/p-) 



-q-(n,T) + g+(r,t), 



(1) 



(2) 



(3) 



where d/dt = d/dt + v ■ V is the Lagrangian derivative, p 
is the mass density, v is the fluid velocity, p is the thermal 
pressure, e = p/(7 — 1) (we use 7 = 5/3, valid for an ideal 
non-relativistic gas) is the internal energy per unit volume, 
g is the acceleration due to gravity, g"*" is the heating rate 
per unit volume, q~ = neniA{T) is the cooling rate per unit 
volume. Tie (rii) is the electron (ion) number density, and A(r) 
is the cooling function. We use a third of the solar metallicity, 
corresponding to a mean molecular weight per particle p = 
p/{nmp) — 0.62 and a mean molecular weight per electron 
Pe = pI (riemp) = 1.18, where n = rie -I- rii is the total number 
density and is the proton mass. 

We model the ICM in spherical (r, 6, </!>) geometry using 
the ZEUS-MP code (Hayes ct al. 2006). We model the dark 
matter halo as a fixed gravitational potential with an NFW 
profile: 



$ = 



IGMq In (1-fr/r,) 
Tb r/rs 



(4) 



where Mo is the characteristic dark matter mass and rs is the 
scale radius (Navarro, Frcnk, & White 1997). The accelera- 
tion due to gravity 



3 



dr 



is directed toward the center of the dark matter halo. For 
simplicity we do not include the contribution of the ICM or 
the galaxies to the gravitational potential. For our fiducial 
cluster model CIO, we take Mo = 3.8 x lO" Mq and = 
390 kpc. We scale the radii according to the self-similar scaling 
for different halo masses; {r, r^} oc Mq^^ . The concentration 
parameter c = r2oo/ra (where r2oo is the virial radius within 
which the mean dark matter density is 200 times the critical 
density of the universe) is fixed to be c = 3.3 for all halo 
masses. With this choice the virial mass is M200 ~ 1.38Mo. 

The equations are solved using operator (and directional) 
splitting into source and transport steps (Stone & Norman 
1992). We combine the heating and cooling rates at each grid 
cell and subcycle the internal energy equation. Cooling for 
simulations without heating (i.e., for cooling-flows with — 
0) is not subcycled but is treated semi-implicitly (see Eq. 7 in 
Paper I); an explicit treatment would lead to computationally 
prohibitive subcycling at almost all timesteps because of short 
cooling times. This less accurate treatment of cooling does 
not affect the results much. All of the simulations are run for 
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5 Gyr, roughly the time between major mergers of halos at 
the current epoch. 



2.1.1 Initial and Boundary Conditions 

We initiahze the ICM in hydrostatic equilibrium with an en- 
tropy profile^ 

K{r) = '^^ Ko + A'loo ( —) " , (5) 

where rioo = 100 kpc for our fiducial run (CIO) and scales 

1/3 

with Mq for other halo masses (Cavagnolo ct al. 2009). 
The entropy profiles that we use cover the range observed in 
clusters. We seed initial density perturbations such that the 
maximum 5p/{p) ~ 0.3, where Sp = p — (p) and (p) is the 
angle-averaged density as a function of r. The density fluctu- 
ations {5p/{p)) are isotropic and homogeneous (in the three- 
dimensional Cartesian sense), and have a one-dimensional 
power spectrum oc k^^^ for 4 ^ kroMt/2% <C 10, where rout 
is the outer radius (200 kpc for the flducial run). Though 
somewhat arbitrary, this choice for the initial perturbations is 
fairly general and enables us to initialize all of our simulations 
identically. For simulations with heating, the resulting quasi- 
steady state does not depend strongly on the initial density 
perturbations. By contrast, the properties of the multi-phase 
gas in cooling-flows do depend on the initial density pertur- 
bations (see Appendix A — this is because TI does not grow 
exponentially in cooling-flows; Balbus & Sokcr 1989). 

We perform our simulations in spherical (r, 9, <j}) coordi- 
nates, where 1 kpc ^ r ^5 200 kpc for the fiducial run CIO, 
^ ^ ^ TT, and 0^0^ 2tt. Again, we adjust the radial 

1/3 

range for difi'erent halo masses such that r oc Mg . We use a 
logarithmic grid in the radial direction and a uniform grid in 
other coordinate directions. We fix the outer electron num- 
ber density to be rie.out ~ 0.0015 cm~^ for all halo masses, 
a value observed for Abell 2199 at 200 kpc (see Johnstone 
et al. 2002; see also Sharma ct al. 2009 for details of the 
setup). The outer radius is close to the scale radius of Abell 
2199, such that the initial temperature increases monotoni- 
cally with radius. We apply infiow (outfiow) boundary condi- 
tions at the outer (inner) radial boundary, such that a steady 
inflow may develop; mass is not allowed to flow out of (in 
to) the computational domain at the outer (inner) boundary. 
Since the cooling time at the outer boundary is much longer 
than 5 Gyr, we do not expect our results to change signifl- 
cantly with a different outer boundary condition.^ We apply 



^ Although an entropy core provides a good fit for most clusters, 
some non-cool-core clusters have a high-density/low-temperature 
'corona' around the central BCG at < 5 kpc (see Sun 2009). We 
do not address these systems in this paper. 

^ The outer boundary conditions do not affect the results for our 
simulations with heating, in which a hot pressure-supported ICM 
is maintained. However, the cooling-flow simulations without heat- 
ing do depend on the outer boundary conditions, as we discuss in 
Appendix A. 



reflective boundary conditions in the 9 direction and periodic 
boundary conditions in the (j) direction. 



2.2 Cooling & Heating 

We use the cooling function A(r) given by Eq. (12) in 
Sharma, Parrish, & Quataert (2010) (corresponding to the 
solid line in their Fig. 1), which is a flt to the equilibrium 
cooling curve of Sutherland & Dopita (1993) for 0.3 solar 
metallicity (see their Fig. 8). Hydrogen starts to recombine 
below ~ 10* K; this abruptly reduces the cooling rate and cre- 
ates a thermally stable "cold" phase. The ICM is locally ther- 
mally unstable at all temperatures above 10* K, but we show 
later (Fig. 3) that the combined effect of cooling, heating, 
and gravity results in a strongly bimodal temperature distri- 
bution in the plasma. We therefore speak of distinct "hot" 
and "cold" phases. 

While the cooling function for the ICM is fairly well- 
known, the heating of the ICM is not yet understood. For 
concreteness we choose a heating rate such that the energy 
deposited per unit volume is constant (g^ oc p"). The simula- 
tions in Paper I demonstrate that the results are not sen- 
sitive to this choice, as long as time is normalized to the 
thermal instability timescale. Note that, while the cooling 
rate {q~ = n^ruK) is a thermodynamic state function, the 
heating rate (r, t) explicitly depends on radius and time. 
Our idealized heating prescription and our feedback heating 
prescription differ in the form of q^{r,t); i.e., the two models 
represent different assumptions about how the feedback en- 
ergy is distributed in the ICM. In our idealized simulations, 
we distribute the heating to instantaneously balance cooling 
at all radii in the plasma and in our feedback simulations we 
make the heating proportional to the mass accretion rate at 
small radii. 



2.2.1 Idealized Heating 

For simulations with idealized heating (Tables 1 & 2) the 
heating rate at any point in the plasma (in erg cm"'' s~^) is 
equal to to the angle-averaged cooling at that radius: 

q+{r,t) = ^r r%-[n{r,t),T{r,t)]smed9d<p. (6) 

This imposes perfect thermal equilibrium at all radii and at 
all times. Although at first Eq. 6 appears overly idealized, 
it provides a useful framework for studying the thermal in- 
stability. Furthermore, the results of Paper I and section 3.5 
show that introducing even large (~ 100% or more) fiuctua- 
tions about the heating rate in Eq. 6 does not significantly 
change the resulting properties of the hot halo (so long as 
global thermal equilibrium is approximately maintained over 
a few cooling timescales). 
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Figure 1. Mass accretion rate (in yr~^) througii tlie inner boundary (at 1 kpc) for different initial entropy profiles as a function of 

time for our cluster runs with idealized heating {left panel) and without heating (i.e., cooling-flows; right panel). In the simulation labels 
the number following 'C is the initial core entropy (Kq in keV cm'^). For simulations with heating the mass accretion rate at late times 
is reduced because the cold filaments are accreted and a hot ICM with txi/tg > 10 is maintained by heating. All of the low entropy 
cooling-fiow runs look similar at late times and reach a cooling catastrophe. The cooling time for runs with Kq > 100 keV cm'^ is long 
enough that a catastrophic cooling-flow does not occur within 5 Gyr even without heating. 



2.2.2 Feedback Heating 

In our feedback simulations we set the heating rate 



47rr3j [1 + (r/rcf)"] 



(7) 



where is a normahzation constant such that the integral 
of the heating rate over the whole volume J q'^-iivr^dr — 
EfbAfinC^, Efb is the feedback heating efficiency, Min is the 
instantaneous mass accretion rate through the inner radius, 
c is the speed of light, is a cutofT radius for heating, 
and a and /3 are parameters determining the heating pro- 
file. We set the cutoff radius equal to the cooling radius 
where the cooling time is 5 Gyr; this gives Vcf = 100 kpc 
for clusters (Mo = 3.8 x 10^" Mq) and ret = 70 kpc for 
groups (Mo = 3.8 x lO" M©). The cutoff radius must be 
larger than the cooling radius to prevent a cooling catastro- 
phe at late times due to cooling beyond the cutoff radius. If 
the cutoff radius is too much larger than the cooling radius, 
however, we require a higher feedback efficiency to balance 
cooling in the core because the feedback deposits a lot of en- 
ergy beyond the cooling radius, where it is not required. The 
parameters a — —1.5 and /? — 3.5 are chosen such that most 
of the energy is deposited near Vd. The choice a — —1.5 en- 
sures that a relatively flat density /entropy core is obtained, 
as observed; for a heating profile steeper (shallower) than this 
we obtain an inner density profile which increases (decreases) 
with radius. The results are not very sensitive to the exact 
choice of the parameters a, /3, and r^. We discuss the depen- 
dence on efb in detail in section 4. 

Since the heating rate in both models (Eqs. 6 & 7) is con- 
stant at all angles for a fixed radius and time {q'^ oc p°), and 
since cooling is sensitively dependent on density (oc p^A[r]), 



fluid elements cooler and denser than average will cool fur- 
ther; though the plasma is globally in thermal equilibrium, it 
is locally thermally unstable. 



2.3 Important Timescales 

In this section we deflne the most important timescales in 
our analysis. The TI growth time (inverse of the exponen- 
tial growth rate) in the isobaric regime is (e.g., see Eq. 8 in 
Sharma, Parrish, & Quatacrt 2010) 



tTi — - 



tc 



3 (2-dlnA/dlnr)' 



where 



^cool 



3 nfcflT 
2 rieniA 



(8) 



(9) 



In Eq. 8 we have assumed volumetric heating {q'^ oc p'') and 
we have ignored thermal conduction which suppresses the in- 
stability at small scales (Field 1965). This approximation is 
justifled because thermal conduction is anisotropic with re- 
spect to the magnetic field (see Paper I for more details). 
For free-free cooling (with A oc T^/^) relevant for clusters, 
tTi ~ (10/9)tcooi. This timescale depends on how the feed- 
back energy is thermalized in the ICM: if instead g"*" oc p, 
the TI timescale is longer, Iti = (10/3)tcooi, because heating 
partially compensates cooling for an overdense blob. We also 
define the free-fall time 



, 2r 
ta - { — 



1/2 



(10) 



where g = d^/dr is the acceleration due to gravity. As we 
discuss later, the ratio tTi/ifi plays a crucial role in the physics 
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of the ICM — if this ratio is small cold filaments condense out 
of the hot phase, resulting in an abrupt rise in the accretion 
rate and in the associated feedback heating. 



3 SIMULATIONS WITH IDEALIZED HEATING 

In this section, we describe the numerical results from simu- 
lations with our idealized heating prescription. These ideal- 
ized heating runs, in which heating at each radius balances 
angle-averaged cooling, are motivated by the observed ther- 
mal equilibrium in cluster cores. The identification that cold 
filaments are due to TI in a globally stable ICM is easily mo- 
tivated in this setup: the ICM evolves only because we seed 
the local TI with density perturbations. Although this ideal- 
ized heating function appears to be finely tuned, we find that 
the results are quite insensitive to large temporal and spa- 
tial perturbations in heating on top of this exact balance (see 
section 3.5 & Paper I). Moreover, the more realistic feedback 
simulations discussed in section 4 reach a quasi-steady state 
similar to these idealized simulations. 

In section 3.1 we discuss the effects of our heating ansatz 
on the mass accretion rate and the cooling to sub-virial tem- 
peratures; we also compare our results to a cooling-fiow. In 
3.2 we discuss our fiducial run CIO in detail. Section 3.3 shows 
that TI and heating leads to a lower limit on fxi/te and the 
core entropy. In section 3.4 we show that TI affects larger 
radii in lower mass halos, and that there is a significant lack 
of plasma below ~ 1/2 — 1/3 of the virial temperature in simu- 
lations with heating. We end this section with a discussion of 
the robustness of our results to large variations in heating on 
top of our thermal balance and to changes in resolution (sec- 
tion 3.5). Table 1 lists the two-dimensional simulations with 
idealized heating (see Eq. 6). Table 2 lists simulations with 
temporal and spatial heating fiuctuations, three-dimensional 
simulations, and simulations at different resolutions. 

3.1 Effect of Heating on Cooling &: Accretion 

To assess the effect of heating on the structure of the cooling 
core we also carried out a few runs without heating [q^ — 0; 
i.e., cooling-flows). In this section we compare the results of 
simulations with idealized heating to cooling-flow simulations. 
The properties of the cooling-flow simulations are discussed 
in more detail in Appendix A. 

3.1.1 Suppression of the Mass Accretion Rate via Heating 

Figure 1 compares the mass accretion rate through the inner 
radius for cluster simulations (Mo = 3.8 x 10^* Mq; pre- 
fixed with C in Table 1) with idealized heating {left panel) 
and without heating {right panel) for different initial entropy 
profiles. The idealized heating runs with initial core entropy 
Kq ^ 30 keV cm^ have a mass accretion rate < 1 Mq yr~^ at 
all times. Runs with lower core entropies (CIO and CI) have 
much larger instantaneous mass accretion rates (up to 100s 
of M0 yr~^), but only for short times; more typically they 



show accretion rates between 1 and 10 Mq yr"'^, in broad 
agreement with accretion rates inferred from X-ray spectral 
modeling (see Peterson & Fabian 2006 for a review on X-ray 
spectroscopy of clusters) and from the star formation rates 
in the central galaxy (e.g., O'Dca et al. 2008). The cooling- 
fiow runs are qualitatively different from the idealized heat- 
ing simulations. After a few cooling times the mass accretion 
rate for the cooling-fiow runs is consistently larger than 100 
Mq yr~^. Only when the central cooling time is longer than 
5 Gyr (as is the case for Kq — 100 keV cm^) does one escape 
this conclusion. Thus, as intended, our idealized heating pre- 
scription prevents a cooling catastrophe even for runs with 
short cooling times. 

Figure 2 shows an alternate view of the fundamental dif- 
ference between cooling-flows and clusters in global thermal 
balance. Here we show the time spent at a given M for sim- 
ulations with and without heating. While the mass accretion 
rate is ~ 100 Mq yr~^ at most times for cooling-flows, simula- 
tions with idealized heating primarily accrete at 1 Mq yr~^. 
Even for the heating runs that show M > 100 Mq yr~^ in- 
stantaneously, the time spent in such a phase is only 10 Myr.'' 
This high accretion rate arises only in low entropy simula- 
tions (see Table 1) and represents cold filaments which con- 
dense out of the hot phase and cross the inner boundary. 
Simulations with intermediate entropies (such as the run C30) 
do not develop multi-phase gas (see Table 1) and never have 
large accretion rates; M <^ 1 Mq yr~^ at all times. Thus, 
large accretion rates in simulations with our idealized heat- 
ing prescription only occurs via cold gas. We show later that 
this result is not sensitive to our choice for the heating. 

The absence of cold filaments in our idealized heating 
simulations with Kg > 30 keV cm'^ appears surprising at first 
because the cooling time in these clusters (« 1 Gyr) is shorter 
than the run-time — since we initialized the ICM with large 
density perturbations, one might expect the TI to become 
strongly nonlinear and to produce multi-phase gas. The cool- 
ing time in these simulations is much longer than the free-fall 
time, however, and gravity significantly alters the saturation 
of the TI in this limit. The density perturbations saturate at 
a small amplitude {Sp/p < 1) because the effective gravity 
felt by slightly overdense blobs {gSp/p) produces strong shear 
which mixes the overdense blobs with the hot phase within a 
cooling time. Thus, overdense fluid elements do not cool all 
the way to the thermally stable phase. By contrast, when the 
TI timescale is shorter than the free-fall time (txi/iff < 1) 
gravity has essentially no effect on the TI. The TI develops 
and saturates over a few cooling times, and only after a free- 
fall time do the cold filaments fall in toward the center (see 
Paper I for a more detailed discussion of the physics of TI in 

^ The Eddington rate for a 10^ Mq black hole is ~ 10 Mq yr"^ (as- 
suming a radiative efficiency of 0.1). Thus, if most of the mass 
accreted in form of cold filaments with M > Impy can reach the 
supermassive black hole, it can power a luminous quasar for a short 
duration (e.g., see Ciotti & Ostriker 1997). Since our simulations 
do not include star formation and accretion physics we cannot ex- 
plore this possibility in detail. 
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Figure 2. The mass accretion rate through the inner radius (e.g., 
in Fig. 1) shows large variations as a function of time, especially 
when cold filaments cross the inner boundary. This plot repre- 
sents the same information as Fig. 1, but shows the time inter- 
val spent accreting with a given mass accretion rate; we use 50 
logarithmically spaced bins in the total mass accretion rate (M) 
from 10~^ to 10^ Mq yr~^. While cooling-flows accrete at large 
rates (~ 100 Mq yr~^) in the cold phase, simulations with heat- 
ing (even the low entropy run CI; not shown in this Figure) spend 
most of their time accreting in the hot phase with M < 1 Mq yr~^. 
The results for our idealized heating model and the feedback heat- 
ing simulations are very similar for both clusters and groups. 
While the mass accretion rate in the hot phase is smaller for 
smaller mass halos, the mass accretion rate is similar whenever 
M > few Mq yr"'^; accretion with M > few Mq yr"-*^ is mainly 
due to the infalling cold gas. 

the presence of gravity). In our setup, cold filaments condense 
out from the hot phase, only if the ratio tTi/tg <^ 10 locally; 
we discuss this more quantitatively below. 

3.1.2 Cooling to Sub-Vmal Temperatures 

The effect of heating is also manifest in the temperature struc- 
ture of the plasma. Figure 3 shows the luminosity per loga- 
rithmic temperature bin for simulations with (CIO) and with- 
out (ClO-cf) heating. We show the simulation with heating 
at two times, one at which there is significant cold gas via 
TI (1-2 Gyr) and one at which there is not (4-5 Gyr). Even 
when filaments are present, the luminosity between ~ 10* K 
and 10'' K is a factor of 10-100 times smaller for the run 
with heating than for a cooling-flow. This is consistent with 
observations that find a striking absence of the low-energy X- 
ray lines expected in a cooling-flow (e.g., Peterson & Fabian 
2006). The luminosity at « 10* K (below which the plasma 




Figure 3. Luminosity per logarithmic temperature 
bin (dL/LogjQT; such that the total luminosity is given by 
the area under the curve) as a function of Log]^QT]j(,v for the 
cooling-flow run (ClO-cf) in the quasi-steady state, for the 
idealized heating run CIO averaged from 1-2 Gyr (during 
which it shows cold filaments), and for CIO averaged from 
4-5 Gyr (when the filaments are absent). We use 1000 bins in 
LogjQTijcVi E^nd the simulation data are averaged in time and fit 
by a high-order polynomial. The results expected from steady 
cooUng-flows (dL/dT = [5/2]fcsTM//imp; e.g., Fabian 1994) 
with M = 250 Mq yr"! and 2.5 Mq yr"'^ (long dashed lines) are 
shown for comparison. 

becomes thermally stable) is also smaller for simulation with 
heating, but only by a factor ~ 3. There is a significant lack 
of plasma at the intermediate temperatures (0.001-1 keV) 
for simulations in thermal balance because gas blobs at inter- 
mediate temperatures are thermally unstable and run away 
either to the hot or the cold phase. Figure 3 also shows that 
the heating simulations do not show any gas below ~ 1 keV 
form 4-5 Gyr (when tTi/tff ^ 10) because there is no cold 
gas produced by TI. 

The luminosity vs. temperature for our coohng-fiow sim- 
ulation in Figure 3 closely matches the analytic cooling-flow 
model between 0.1 to 1 keV {dL/dT = [5/2]kBTM/fimp; 
e.g., Fabian 1994), with a cooling rate M = 250 Mq yr~^ (a 
value similar to the accretion rate in cooling-flow simulations; 
see Table 1). Without a heat source, the plasma cools al- 
most monolithically; our cooling-flow simulations thus possess 
roughly equal mass per cooling time interval. When heating 
balances cooling, however, the ICM does not cool monolithi- 
cally. Slightly overdense cool blobs are therefore able to mix 
with the hotter ambient plasma as they fall in. Thus cooling 
from 1 to 0.1 keV follows a random walk in cooling/heating 
rather than cooling monolithically (Soker & David 2003 en- 
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Figure 4. Snapshots of the electron number density in the inner 
10 kpc for the fiducial simulation (CIO) at different times around 
when cold filaments condense out of the hot phase. The limits on 
the electron density contour plots only go from 0.01 to 1 cm~^ 
but the coolest/densest plasma is extremely cold (~ 10* K) and 
dense (~ 100 cm~^). Cold filaments condense out of the hot phase 
and eventually fall in toward the center; however, the cold gas does 
not just free-fall radially inwards and can even move outwards in 
radius for some time. 



visage a similar mechanism of a 'moderate cooling-flow') and 
the slope of the luminosity function is much steeper. More- 
over, the amount of gas cooling through lower X-ray tempera- 
tures is smaller, as observed, because of reheating. Therefore, 
in simulations with heating the net cooling rate below 1 keV, 
which is approximately equal to the average mass accretion 
rate in the cold phase at small radii, is 10-100 times smaller 
than a cooling-flow.'' However, even in the presence of heat- 
ing, once gas manages to cool below ^ 1 keV, it runs away 
all the way to 10* K. Thus our idealized heating run resem- 
bles a much-reduced cooling-flow below ~ 1 keV (cooling is 
reduced by ~ 100; compare with the cooling-flow prediction 
for M = 2.5 Mq yr-i). 

The cooling- flow prediction using M = 250 M0 yr"'^ 
in Figure 3 matches our cooling only simulation for tem- 
peratures above 0.1 keV. However, the results deviate at 
lower temperatures for two reasons: first, a cooling-flow with 
large density perturbations, as is the case in our simulations, 
can give rise to blobs at ~ lO** K cooling out of the soft 
X-ray emitting gas, whereas a steady cooling-flow is much 



* The exact suppression of luminosity relative to a cooling-flow 
depends on the ratio tTi/tff ■ For a small tTl/%) E.g., at early times 
in the low entropy run CI the luminosity below 1 keV is closer 
to the cooling- flow value. For larger txi/iff, the luminosity below 
1 keV can be vanishingly small, as shown by the CIO results at 
4—5 Gyr in Figure 3. 



smoother (see Appendix A for more details); second, our reso- 
lution is not high enough to obtain robust statistics for plasma 
at ~ 10^ — 10^ K. A plasma blob cooling isobarically from 
lO'^ K shrinks by a factor ~ 10* in volume by the time it 
reaches the stable phase.'' Thus cold filaments are not spa- 
tially well-resolved even in our highest resolution simulations! 
However, we are able to get reasonable statistics by combining 
simulation data over a long time interval. Very high resolution 
simulations of the core to quantitatively calculate emission in 
different temperature bins will be carried out in the future. 



1.5 3.2 Cool-Core Clusters in Thermal Balance 

Figure 4 shows snapshots of the inner radii in our fiducial 
idealized heating run (CIO) at times when cold filaments con- 
2 dense out of the hot phase. The top-left panel shows a snap- 
shot right before the condensation of a cold 'filament.' Once 
the over-dense blob cools below ~ 1 keV, the cooling rate 
becomes faster than the turbulent mixing rate and the blob 
quickly cools to the thermally stable temperature. This newly 
born cold filament has a non-zero velocity because of the mo- 
tion of the overdense blob from which it formed; cold filaments 
can thus have instantaneous velocities directed away from 
the center of the potential. Further cooling of the hot phase 
close to the cold filament seeds condensation of secondary fil- 
aments. Since filaments are born with non-zero velocity, they 
move through the ICM core before leaving through the inner 
boundary after ~ few free-fall times. The cold blobs at the 
thermally stable temperature (~ 10* K) are in rough pres- 
sure equilibrium with the hot phase. The dramatic difference 
in density (and temperature) of the cold and hot phases, and 
the lack of plasma at intermediate densities (temperatures) 
clearly points to local TI as the origin of these cold filaments. 
The production of cold filaments ceases once the density of 
the hot phase is depleted such that tTi/tg > 10; the bottom- 
right panel of Figure 4 shows the core in a quasi-steady state 
without filaments. 

Multiphase structure also develops in cooling-flows with 
large density perturbations because of a faster cooling rate 
in overdense blobs (see Appendix A for a detailed analysis). 
Whereas the formation of multi-phase gas in our simulations 
with idealized heating depends only on the ratio tTi/ta, pro- 
duction of multi-phase gas in cooling-flows depends also on 
the amplitude of the density perturbations. 

Figure 5 shows the time- (over 1 Gyr intervals) and 
angle-averaged profiles of entropy (left panel) and the ratio 
iTi/tff (right panel) for the fiducial run from to 5 Gyr. The 
cold filaments of the kind seen in Figure 4 show up as troughs 
in the entropy and txi/tff profiles at 1-2 and 2-3 Gyr. Hav- 
ing condensed from the hot phase, these dense filaments fall 
through the inner boundary at roughly the free-fall rate. By 
3 Gyr the hot ICM has both been heated and lost sufficient 



^ For a steadily cooling blob M/t^^ol ~ constant, irrespective of 
the temperature bin, where M is the mass in a given temperature 
bin. Thus, the volume occupied by a cooling blob in an isobaric 
cooling-flow is reduced by ~ T^/A. 
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Figure 5. Time- and angle- averaged entropy {left; {rieT) / (ue)^^^ , where (} represents averaging in and (/>) in keV cm'^ and iTi/tff (fight) 
as a function of radius for the fiducial run CIO. Solid black line represents the initial profile. The other lines are the profiles averaged over 
a 1 Gyr timescale at different times. Sharp low entropy (and low txi/tff ) troughs seen from f— 2 Gyr and 2-3 Gyr are because of infalling 
cold gas filaments, which are also responsible for the large accretion rates seen in Figures 1 & 6 during the same time intervals. 



mass due to "drop-out" of cold filaments that tTi/ts > 10 
everywhere. There are no more cold filaments produced once 
ixi/iff > 10, even though the cooling time in the core is less 
than 1 Gyr. This is also reflected in the entropy proflle as a 
rise in the central entropy to Kq k. 20 keV cm'^ in the final 
state. As mentioned already, cold filaments condense out of 
the hot phase only when the ratio tTi/iff <~ 10, or equivalently 
when the core entropy Kq <: 20 keV cm^. 

The cold filaments falling toward the center cause a sud- 
den increase in the mass accretion rate at small radii. Fig- 
ure 6 shows the mass accretion rate through the inner ra- 
dius (1 kpc) in different phases for the fiducial run: cold 
phase with T < 0.01 keV, intermediate phase with 0.01 keV 
< T < 0.1 keV, and hot phase with T > 0.1 keV. These 
temperature ranges are chosen somewhat arbitrarily, but the 
results do not depend sensitively on the exact temperature 
values because there is very little plasma between the stable 
phase (at ~ 10* K) and about one-third of the virial tempera- 
ture (see section 4.3). The big spikes in M for the fiducial run 
CIO in Figure 1 coincide with cold filaments (seen in Fig. 4) 
falling through the inner boundary. The mass accretion rate 
in the hot and intermediate phases is < 10 M© yr"'^ at all 
times. However, the mass accretion rate in the cold phase can 
reach up to 1000 yr~^ for short times due to cold dense 
filaments crossing the inner boundary. The infalling cold fila- 
ments have a sheath of plasma at intermediate temperatures 
which contributes to the small mass flux at intermediate tem- 
peratures during the infall of cold filaments. 



3.3 Self-Regulation of the ICM core 

Having discussed how the timescale ratio tTi/tff influences 
the accretion rate and the amount of cold gas in the ICM, 
we now discuss how these effects in turn influence the ratio 



r ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' '1 




time (Gyr) 

Figure 6. Mass accretion rate through the inner boundary (1 kpc) 
for the fiducial run (CIO) as a function of time. The accretion 
rate in different temperature phases is indicated by different col- 
ors/styles. The instantaneous accretion rate in the cold (T < 
0.01 keV) phase can be large, reaching > 100 Mq yr~^, but lasts 
only for short bursts (see Fig. 2). Most of the time accretion occurs 
in the hot (T >0.1 keV) phase. Accretion of plasma at intermediate 
temperatures (0.01 keV < T < 0.1 keV) coincides with the infalling 
cold filaments, but is subdominant. The green short-dashed line in- 
dicating accretion at intermediate temperatures is obscured by the 
blue long-dashed line showing accretion in the cold phase. 
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iTi/tft in the cool-core of the ICM. We show in this section 
that the local TI tends to maintain clusters close to the critical 
threshold for forming filaments such that txi/tft ^ 10. The 
condition for forming cold filaments (iTi/tff < 10) can be 
expressed in terms of entropy and the free-fall time as 

2/3 



K <20 keV cm^ 



tfi 



^SOMyryj ' ^^^^ 

where we have assumed a heating rate which is constant per 
unit volume {q~^ oc p''). Thus our thermal instability model 
can explain the observational lack of cold gas indicators in 
clusters with Kq > 30 keV cm^ and the observed over- 
abundance of cool-core clusters at Kq ~ 10 keV cm^ (e.g., 
Fig. 1 in Cavagnolo et al. 2008). 

Figure 7 shows the entropy and tTi/tfj profiles at dif- 
ferent times for a cluster with an initially low entropy core 
CI {Ko — 1 keV cm ). Cold filaments quickly condense out of 
the hot phase since the ratio tTi/tg < 10 everywhere within 
20 kpc. The infalling cold filaments leave the inner boundary 
after roughly a free-fall time. The entropy of the remaining 
plasma, maintained in hydrostatic equilibrium by our ideal- 
ized heating prescription, thus increases. For extremely low- 
entropy initial conditions the central entropy increases very 
quickly because of a short cooling time; this happens in less 
than 1 Gyr for the run CI (see Fig. 8). After a few TI times, 
the density profile (and the entropy profile as well) adjusts 
such that the ratio tTi/ta is close to the critical value of fs 10. 
This adjustment of the inner ICM because of TI and infalling 
cold filaments plays a very important role in determining the 
density profile of the hot plasma in halos of different masses, 
as we discuss further in section .3.4. 

The average entropy at small radii from 4-5 Gyr for 
the run CI is slightly lower than the inner entropy at late 
times for the fiducial run (see Fig. 5). This slightly lower en- 
tropy (and txi/ifi smaller than the critical value) results in 
another episode of filaments forming and falling through the 
inner boundary at 4 Gyr, as seen in the left panel of Figure 1. 



Thus, tTi/ts < 10 is a robust criterion for the presence of 
filaments, irrespective of the initial conditions. Moreover, the 
critical ratio tTi/ta is self-regulated to be ^ 5 — 10 because of 
heating and mass drop-out in the form of cold filaments. Also 
notice that the inner core at late times does not have a fiat 
entropy profile, though the profile is shallower compared to 
the larger radii. While the exact entropy profile at small radii 
depends on the details of heating and cooling, an isentropic 
core appears to be a good approximation because the cool 
low-entropy plasma blobs leave from the inner boundary and 
the hot plasma is mixed throughout the core. 

Figure 8 shows the minimum of the ratio tTi/te as a 
function of time for the cluster runs CI, CIO, and C30 (and 
the group run G50 discussed in section 3.4). As can be seen in 
Figures 5 & 7, the minimum in tTi/ta increases due to mass 
drop-out in the form of cold filaments for runs CI and CIO. 
Whenever tTi/tg <^ 10 in Figure 8, cold filaments condense 
out of the hot phase and fall in through the inner bound- 
ary, and are seen as spikes in the mass accretion rate in 
the left panel of Figure 1. The minimum of the ratio tTi/tg 
is larger than the critical value (10) for the run C30 (with 
Ko — 30 keV cm ) and there are no cold filaments (e.g., see 
Table 1). Figure 8 shows that this ratio (and the core entropy 
Kq) does not evolve much with time, even though the cooling 
time is ~ 1 Gyr. The mass accretion rate for this run is much 
smaller than for the cooler clusters because massive cold fil- 
aments cannot condense out of the hot phase (see left panel 
of Fig. 1). Although the cooling time is less than the cluster 
age for run C30 there is no cooling-fiow because of our heat- 
ing prescription (compare C30 and C30-cf in Fig. 1). Thus, 
even clusters which are unlikely to show cold filaments (since 
tTi/tft > 10), but which have cooling times shorter than the 
cluster age, require heating roughly balancing cooling to pre- 
vent a massive cooling-flow. 

The rapid adjustment of cool-core clusters to tTi/tg > 10 
via mass drop-out and heating, and the suppression of M (and 
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Figure 8. The minimum of the ratio tTl/% as a function of time 
for various runs (evaluated from angle- averaged radial profiles of 
plasma with temperature > 0.1 keV sampled every 50 Myr). Cold 
filaments condense out and result in large accretion rates whenever 
the ratio fxi/tff < 10 (compare with the left panel of Fig. 1). This 
in turn leads to significant heating, driving tTi/tff to be > 10. 



hence feedback) if tTi/tff > 10 suggest that cool-core clus- 
ters will spend most of their time close to the threshold 
for forming cold filaments, so long as heating is related to 
the amount of cold gas produced by TI (as is the case for 
AGN/supernova feedback). This physics can quantitatively 
account for the observed excess of clusters with central en- 
tropies around 10 keV cm'^ (e.g., Cavagnolo ct al. 2008; 
Pratt ct al. 2010). In addition to cool-core clusters affected 
by cooling and feedback, there are also non-cool-core clusters 
which have high entropies because of strong feedback heating 
and/or mergers. These non-cool-core clusters, relatively un- 
affected by cooling, maintain their high entropies and form a 
class distinct from the cool-core clusters. 

3.4 Comparison of Different Halo Masses 

We have also carried out simulations with our idealized heat- 
ing prescription in lower mass, group-sized halos (Mo = 
3.8 X lO" M0), in massive groups with Mo = 10^'' Mq, and in 
more massive cluster halos (Mo = 10^^) — these appear in Ta- 
ble 1 with the prefixes 'G,' 'MG,' and 'MC,' respectively. The 
number after the prefix represents the entropy of the equiv- 
alent fiducial-mass halo, assuming self-similar scaling. Thus, 
runs G300 and C300 have identical density profiles, but the 
temperature in run G300 is lower by a factor of 10^''^ and the 
core entropy in G300 is Kq = 300/10^^^ = 64.6 keV cm^ 

At a given scaled entropy, cooling becomes stronger in 



lower mass halos because of the corresponding decrease in 
the virial temperature (fcooi cx: T^^"^ jn oc M^''^). On the other 
hand, the free-fall time (iff oc [r^/Monc]^''^) at a given (scaled) 
radius is roughly independent of the halo mass because dark 
matter, which dominates gravity, is self-similar. Consequently, 
groups develop multi-phase gas at a higher scaled entropy 
than do clusters; e.g., while cold filaments do not condense 
from the hot phase in the cluster run C30, filaments are seen 
in the group run G50. 

The group run G50 is, in fact, initially very close to the 
threshold for forming multi-phase gas. Figure 8 shows the 
minimum of the ratio tTi/ift as a function of time for this 
simulation. It does not show cold filaments from 0-1 Gyr, but 
because of slightly low entropy material falling in from larger 
radii, it crosses the threshold for forming cold filaments after 
1 Gyr. The group also undergoes a second, stronger episode of 
thermal instability and again forms filaments after « 3 Gyr; 
this depletes the core of much of its gas and raises the entropy 
in the hot phase much beyond the critical threshold. 

Figure 2 shows that the mass accretion rate in the hot 
phase is smaller in groups than in clusters by a factor of ~ 
10 (their mass ratio; see also Table 1), but that the mass 
accretion rate in the cold phase (which accounts for any M ^ 
few Mq yr~^) is nearly independent of the halo mass. Thus, 
the amount of gas which reaches the cold phase is similar in 
our group- and cluster-mass halos and these halos therefore 
receive similar amounts of feedback energy. The mass of the 
hot phase is lower in groups than in clusters, however, and the 
effect of this feedback is therefore much stronger. Thus, even 
in our idealized heating model, heating of the hot phase can 
overwhelm cooling in lower mass halos, resulting in a much 
lower density (and larger tTi/tff) than the critical threshold 
for filament formation. The last 2 Gyr in Figure 8 illustrate 
this overshoot. 

3.J^.1 Bigger Cool-Cores for Lower Mass Halos 

Figure 9 shows the angle averaged entropy and electron num- 
ber density in the hot phase (T > 0.1 keV; time averaged from 
4-5 Gyr) as a function of radius (scaled to the virial radius) 
for simulations with different halo masses; only halos prone 
to the formation of cold filaments according to our ixi/iff cri- 
terion (runs MCI, CIO, MGIO, & G50) are shown. The ICM 
with tTi/iff <, 10 develops cold filaments which remove mass, 
and bring Iti /tg , entropy, and density in the core close to the 
critical value. Thus, the core density (entropy) for each halo 
in Figure 9 represents a limit below (above) which cold fila- 
ments cannot condense out for a given halo mass. While the 
temperature decreases with a decreasing halo mass, the core 
radius (where the ratio tTi/tft oc T^^"^ /{nta) is the smallest) 
increases, assuming a similar density profile outside the core 
radius. Thus the free-fall time tg oc (r/r2oo)'^^^ at the core 
radius increases with a decreasing halo mass. These trends in 
temperature and tg at the core radius roughly cancel, lead- 
ing to an entropy threshold for filament formation (Eq. 11) 
that is similar for groups and clusters. This can be verified 
from Figure 9 which shows a relatively constant entropy for 
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Figure 9. The angle- and time-averaged (4-5 Gyr) electron num- 
ber density (rie incm~^) and entropy (keV cm^) as a function 
of radius (scaled by the virial radius) for the idealized heating 
runs (MCI, CIO, MGIO, & G50; solid lines) and for the runs with 
feedback heating (C10-fb-2em3 with efb = 2 X IQ-^ &; G50-fb-7em5 
with Efb = 7 X 10~^; dotted lines). Only the hot X-ray emitting 
plasma with T > 0.1 keV is used to calculate the average radial 
profiles. Both the core density and temperature decrease with a 
decreasing halo mass; therefore the core entropy is relatively inde- 
pendent of the halo mass. 



different halo masses; as a result, the core entropy for groups 

2/3 

is larger than the self-similar prediction K oc Mq , as is ob- 
served (e.g., Ponman, Caimon, & Navarro 1999; Sun ct al. 
2009). 



3.4.2 Lack of Gas below To/3 

In section 3.1.2 we showed that for our cluster model in ther- 
mal balance only a tiny fraction of the plasma cools below 
« 1 keV. Here we discuss the cooling of the hot gas in different 
halo masses. Figure 10 shows the luminosity (dL/d log jg T) as 
a function of temperature for simulations with different halo 
masses which show cool filaments (MCI, CIO, MGIO, & G50; 
these are all 'cool-core' halos). The axes are scaled accord- 
ing to the self-similar prediction. The profiles are obtained 
by averaging over 1 Gyr coinciding with the last episode of 
condensation due to TI. For all halos the luminosity drops 
quite abruptly below the temperature corresponding to the 
peak luminosity. Moreover, the shape of dL / dLog-^^fT vs. tem- 
perature below the peak temperature is quite similar, and is 
consistent with the observational fact that there is negligible 
luminosity below 1/2 — 1/3 of the ambient cluster tempera- 
ture (e.g., Peterson ct al. 2003; Sanders et al. 2010). The line 
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Figure 10. Luminosity per logarithmic temperature 
bin (dL/Logj^pT) as a function of temperature for different 
halo masses; the data are averaged over a 1 Gyr interval when 
the models have cold filaments. The axes are scaled according to 
the self-similar prediction for luminosity (oc Mg^^) and tempera- 
ture (oc Mq^^). The normalization of both axes is relative to the 
cluster mass halo (Mc = 3.8 X 10^* Mq). 



corresponding to the group G50 in Figure 10 shows a small 
bump at double the ambient temperature because some of the 
gas is heated above the virial temperature; this is consistent 
with the late-time rise in the ratio txi/iff for G50 in Figure 8. 



3.4.3 The Required Heating Efficiency 

The last columns of Tables 1 & 2 list the required average 
heating efficiency e — Q'^ /Mc? corresponding to our heating 
prescription, where Q"*" = J^"°°^ (r)A7vr^dr {q~^ is given by 
Eq. 6), rin is the inner radius, rcooi is the cooling radius where 
the angle-averaged cooling time is « 5 Gyr, and M is the av- 
erage accretion rate through the inner radius. In calculating 
the efficiency we do not integrate over the whole domain, al- 
though we do apply our heating prescription at all radii, be- 
cause the radii where the cooling time is longer than 5 Gyr do 
not require heating to maintain the pressure support. Thus, 
the efficiency is zero for the high entropy runs with cooling 
times longer than 5 Gyr. Among cluster runs (prefixed by C), 
the required efficiency (e) peaks at ~ 5 x 10"'' for the low- 
est core entropy run which does not show cold filaments (run 
C30). Even for other halo masses, the required heating effi- 
ciency roughly peaks for the lowest entropy run not showing 
filaments. The required efficiency is smaller once cold fila- 
ments can condense because while the global cooling rate is 
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similar to the higher entropy runs with no filaments, the mass 
accretion rate is dramatically larger. 

The required heating efficiency for runs at the threshold 
for forming cold filaments (MCI, CIO, MGIO, and G50) de- 
creases with a decreasing halo mass; e.g., while e for the run 
CIO is 10~^, it is 7.4 x 10"^ for G50, more than ten times 
smaller. This is because the mass in the hot phase decreases 
by ~ 10 from CIO to G50, but the mass accretion rate (pri- 
marily in the cold phase) is about the same (Fig. 2; see also 
Table 1). Thus, a smaller feedback efficiency is sufficient to 
balance cooling for a smaller mass halo. For a realistic feed- 
back scenario this means that while AGN feedback is required 
for group and cluster mass halos (supernova feedback is quite 
inefficient, e < 10~^), supernova feedback may be sufficient 
for Galactic scale halos that have an even smaller mass in the 
hot phase. A similar conclusion is reached for our simulations 
with feedback heating discussed in section 4. One implication 
of this result is that, if AGN heating with a fixed feedback 
efficiency is responsible for heating cool cores, the lower mass 
halos can be heated much more than more massive ones. 



3.5 Robustness of the Idealized Heating Ansatz 

3.5.1 Spatial & Temporal Perturbations on Heating 

To test the robustness of our results we also carried out sim- 
ulations with spatial and temporal perturbations in heating 
on top of thermal balance at every radius. Simulations with 
rather large variations about this balance are qualitatively 
similar to the idealized simulations (Table 2; see also Pa- 
per I). The run ClO-plO is similar to the fiducial case but 
with temporal perturbations in the heating rate of a factor 
of ~ 10 at each grid cell (on top of given by Eq. 6). The 
perturbations, chosen from a uniform distribution centered 
at zero, have a correlation time of 100 Myr (the results do 
not depend sensitively on the correlation time). The entropy 
and tTi/iff profiles are similar to those for the fiducial run in 
Figure -5. Table 2 also shows that the higher entropy run with 
perturbations C30-pl0 does not show cold filaments, consis- 
tent with our idealized run C30. Thus, our picture of TI holds 
even with large perturbations in the heating rate because on 
average the plasma is still in thermal equilibrium. 

We also carried out simulations (ClO-nzl, C10-nz4, etc.) 
in which the balance of heating and cooling is imposed on 
larger scales, rather than at every radius as is the case in the 
idealized simulations we have focused on; e.g., in run ClO- 
nzl6 the heating rate in each of 16 logarithmically spaced ra- 
dial zones is calculated by averaging cooling over each zone. 
From these simulations we find that if heating balances cool- 
ing on scales smaller than the cooling scale height (the scale 
over which the cooling rate changes significantly) then the 
results are similar to the idealized heating runs. We thus con- 
clude that, as argued in Paper I, many of the details of heat- 
ing are not important for the thermal evolution of the halo, 
as long as the heating stabilizes the ICM against a cooling 
catastrophe. 



3.5.2 Three-Dimensional Runs & Convergence Studies 

Table 2 lists three-dimensional (256 x 128 x 32) simulations 
of clusters CIO and C30. The two- and three-dimensional 
simulations agree in that the runs C30 and C30-3D with 
iTi/iff > 10 do not show cold filaments but the runs with the 
ratio tTi/tft < 10 (CIO and C10-3D) show filaments. The mass 
in the cold phase, mass accretion rate, etc. are not identical in 
two and three dimensions because there are non-axisymmetric 
initial perturbations in the three-dimensional simulations. 

Table 2 also lists two-dimensional simulations but with 
higher and lower resolution than the fiducial case. The lower 
core entropy (CIO, ClO-hlf, ClO-dbl) cases show cold fila- 
ments for all resolutions. Similarly, the higher core entropy 
runs (C30, C30-hlf, C30-dbl) do not show cold filaments for 
any resolution. The threshold for forming multi-phase gas and 
the integrated quantities such as the average mass accretion 
rate are roughly independent of the resolution. We showed in 
Sharma, Parrish, & Quatacrt (2010) that thermal conduc- 
tion must be included in order to achieve strict convergence. 
Since we cannot resolve the Field length in the cold phase 
with a realistic cooling function using present computational 
resources, we necessarily sacrifice strict convergence. 



4 SIMULATIONS WITH FEEDBACK HEATING 

In this section we describe the feedback simulations in which 
the integrated (isotropic) heating rate is proportional to the 
instantaneous mass accretion rate at the innermost radius. 
This feedback model is clearly only a rough approximation 
to heating by a central AGN, but as we shall show, it pro- 
vides a useful demonstration that many of the results in the 
previous section are robust and do not depend sensitively on 
the precise details of the heating. Although the exact mecha- 
nisms which heat the ICM are not fully understood, possible 
scenarios involve jets and bubbles launched by the central 
AGN. These processes are ultimately fed by accretion, how- 
ever, and their power is likely to be proportional to the mass 
accretion rate. Feedback heating by jets, bubbles, or star for- 
mation likely lags the instantaneous accretion rate by up to 
~ 100 Myr, the dynamical time at the core radius. Moreover, 
the feedback efficiency is likely to be time variable and may 
depend on the mass accretion rate. More realistic modeling 
of accretion and launching of jets is outside the scope of this 
paper and will be investigated in the future. 

Our feedback simulations are initialized far from equi- 
librium, with the heating rate = but the cooling rate 
q~ > 0. Nonetheless, after a transient phase, they are able 
to reach an approximate thermal equilibrium where average 
heating in the core roughly balances average cooling. Thus, 
the results from our idealized heating simulations carry over 
to these more realistic feedback simulations. Table 3 lists the 
runs with feedback heating using g"*" in Eq. 7. In section 4.1 we 
describe the feedback simulations for clusters. In section 4.2 
we show that feedback simulations for groups require lower 
feedback efficiency to achieve thermal balance. Section 4.3 
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time (Gyr) 

Figure 11. Total mass accretion rate (in Mq yr~^) through the 
inner radius (1 kpc) for the cluster runs with feedback heating. 
Different feedback efficiencies are labeled. Large spikes in the ac- 
cretion rate correspond to infalling cold filaments crossing the inner 
boundary. The duty cycle for cold filaments and AGN feedback is 
shorter for a smaller feedback efficiency. 

compares the results from our idealized heating simulations 
discussed in the previous section with the feedback runs. 

4.1 Feedback Heating in Clusters 

Figure 11 shows the mass accretion rate as a function of time 
for the cluster simulation CIO with our feedback heating pre- 
scription. We show simulations with feedback efficiencies eft 
= 0.007, 0.004, 0.002, 7x 10"" (see Eq. 7). As discussed above, 
all simulations are initially out of thermal equilibrium and the 
accretion rates spike around 0.2 Gyr, or roughly the cooling 
time in the core. If the feedback efficiency is large enough (e.g., 
Efb = 0.007), the ICM approaches a steady-state with a low 
M and a low core density. The ratio of the crucial time-scales 
tTi/ta is greater than 10 at all radii in this state, and cold 
filaments therefore do not form at later times. If the efficiency 
is too low (e.g., efb = 7 X 10~*), however, a larger mass flux 
is required to maintain thermal equilibrium in the core; such 
high accretion rates can only be supplied by cold filaments. 
The evolution of these simulations is therefore sporadic, punc- 
tuated by episodes of strong cooling or heating. For interme- 
diate feedback efficiencies (efb = 0.002, 0.004), the ratio of 
time-scales tTi/tg is typically larger than 10, but occasion- 
ally drifts to lower values. Filaments in these simulations are 
a rare, transient phenomenon. 

In all cases, cold gas condenses from the ICM when (and 
only when) the ratio tTi/tg < 10. These infalling cold fila- 



ments drive strong feedback and can overheat the ICM, as 
we discuss below. The feedback efficiency sets the duty cycle 
for the production of the multi-phase gas and for feedback 
heating (the duty cycle is shorter for a smaller feedback effi- 
ciency), but not the critical threshold for the TI. The excur- 
sions about the critical threshold can be large, however, and 
depend on the feedback efficiency. The results are not sensi- 
tive to the initial core entropy and the initial transients are 
wiped out on the core cooling timescale. 

Figure 12 shows profiles of entropy and of the ratio 
tTi/tft for a simulation with an intermediate efficiency efb = 
0.002 (C10-fb-2em3). Cooling initially dominates heating and 
causes the ratio txi/ifi to drop below the critical threshold of 
~ 10; the TI thus creates cold filaments which reach the 
center around 0.2 Gyr and lead to feedback heating. Heating 
continues until tTi/tg > 10 everywhere and, by 1-2 Gyr, the 
core entropy has significantly increased. With this intermedi- 
ate efficiency, however, the hot-phase accretion is too small 
for the feedback to completely balance cooling. The core cools 
slowly but secularly and again crosses the tTi/tg <. 10 thresh- 
old after another 1-2 Gyr. This leads to a second episode of 
accretion and heating. This pattern — with the ICM main- 
taining rough thermal balance and producing multiphase gas 
whenever Iti < 5 — 10 iff — would likely continue if we ran the 
simulation further. 

The bump around 60 kpc in the profiles in Figure 12 
result from the particular choice of the heating parameters in 
Eq. 7. Our aim here is not to produce a profile which perfectly 
matches observations by varying all the parameters in the 
heating function, but rather to obtain something reasonable 
without much fine tuning. 

4.2 Feedback Heating in Groups 

Figure l.'"i shows the mass accretion rate through the inner 
radius as a function of time for group mass halos (Mq = 
3.8 X lO^'' Mq) with different feedback efficiencies. The cutoff 
radius (rcf) in Eq. 7 is chosen to be 70 kpc. Other param- 
eters in the heating function (Eq. 7) are chosen to be the 
same as before {a — —1.5 and /3 — 3.5). The qualitative be- 
havior of M as a function of time is similar to the clusters in 
Figure 11 in that accretion in the cold phase becomes progres- 
sively more important as the feedback efficiency is reduced. 
Notice, however, that the transition from a large accretion 
rate in the cold phase to a state with accretion primarily in 
the hot phase takes place at a lower feedback efficiency. For 
clusters this transition happens around efb ~ 0.002 but in 
groups it occurs for efb « 7 x 10~^. This conclusion was also 
reached for simulations with idealized heating and has im- 
plications for possible sources of feedback in groups/clusters 
and galactic halos. 

Figure 14 shows the entropy and tTi/tg profiles for the 
group run with efb = 7 X 10"^ (G50-flD-7em5). The qualita- 
tive behavior is very similar to the cluster simulations with 
intermediate feedback. There is an initial spike in accretion 
because heating is small at early times. After an early episode 
of filament condensation and heating the critical ratio iTi/tff 
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Figure 12. Time- and angle-averaged entropy {left; (ueT) / {rie)^^^) in keV cm^ and iTi/tff (Tight) as a function of radius for the cluster 
run with feedback heating (eft = 0.002; run C10-fb-2em3 in Table 3). Solid black line represents the initial profile. The other lines represent 
the profiles averaged over a 1 Gyr timescale at different times. 




lime (Gyr) 

Figure 13. Total mass accretion rate (in Mq yr^^) through the 
inner radius for group models (G50) with feedback heating. Dif- 
ferent feedback efficiencies are labeled. Large spikes in accretion 
correspond to infalling cold filaments crossing the inner boundary. 



4.3 Comparison to Idealized Heating Simulations 

Figure 15 illustrates the feedback regulation process under 
each of our heating prescriptions. Here we plot the minimum 
of the angle-averaged tTi/ts as a function of time, taking 
only the X-ray emitting plasma into account; this quantity 
determines whether the ICM develops multi-phase structure, 
and its temporal evolution depends qualitatively on the heat- 
ing mechanism. With our 'idealized' heating prescription, the 
heating responds both locally and instantaneously to cooling 
and it keeps the ICM very close to the threshold for filament 
formation, tTi/tft ~ 10. Our 'feedback' heating prescription is 
not so tailored and it, by contrast, only balances cooling in an 
average sense. Consequently, the departures from the critical 
value of tTi/tg are much larger in our feedback simulations 
than in our idealized simulations; this leads to the episodic 
condensation of cold filaments discussed in section 4.1. 

Though simulations with our idealized heating prescrip- 
tion are always in global thermal equilibrium, our feedback 
simulations are initialized very far from equilibrium. As dis- 
cussed above (section 4.1), this leads to a period of rapid, 
cold- mode accretion (around 0.2 Gyr in Figs. 11 & 13), fol- 
lowed by a strong feedback event.'' This is seen as a sudden 
increase in min(tTi/tff), after which heating is suppressed be- 
cause mass accretion occurs only through the much slower 
hot-mode (see Figs. 11 and 13). As the core slowly cools, the 
ratio tTi/ts again falls below the threshold of ~ 10, lead- 
ing again to a period of cold-mode accretion and a strong 



becomes > 10 and there is a period of ~ 1.5 Gyr (from 1.5 
to 3 Gyr) without multi-phase gas. There is a second episode 
of filaments at ~ 3.5 Gyr. One important difference from the 
cluster runs is that the relative size of the core affected by 
cooling and feedback is much larger for groups than for clus- 
ters (see also Fig. 9). 



° The delay between condensation of cold filaments and large en- 
ergy injection due to feedback is of order the free-fall time of fil- 
aments from the radius at which they form. Also notice that the 
secular cooling time in the state without filaments is longer for 
groups than for clusters because of a longer cooling/free- fall time 
in a bigger core (see Fig. 9). 
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Figure 14. Time- and angle-averaged entropy (left; {ueT) / (rie)^^^) in keV cm^ and iTi/tff (right) as a function of radius for the group 
run with feedback heating (efj, = 7 X 10~^; run G50-fb-7em5 in Table 3). Solid black line represents the initial profile. The other lines are 
the profiles averaged over a 1 Gyr timescale at different times. Cooling and feedback affect a larger radius in lower mass halos because of 
the shorter cooling time (e.g., compare this Figure with Figs. -5 & f2). 



feedback event. Eventually the feedback simulations settle in 
rough thermal balance similar to the idealized heating runs. 

There are several additional similarities between our ide- 
alized heating simulations and the feedback simulations. The 
luminosity at sub-virial temperatures (as in Fig. 3) for feed- 
back simulations looks very similar to the idealized heating 
simulations. In both cases cooler gas exists only when the 
ratio tTi/tff < 10, and cooling to sub-virial temperatures is 
suppressed by a factor of ~ 100 as compared to the cooling- 
flow simulations, in agreement with observations of cool-core 
clusters (e.g., Peterson & Fabian 2006). Figure 2 shows that 
the time spent accreting at a given M is very similar for our 
idealized heating simulations and our feedback simulations, 
for both groups and clusters. The fact that the more realistic 
feedback heating runs are much closer to our idealized heat- 
ing runs than to a cooling flow and consistency with several 
cluster observations strongly suggest that our models with 
idealized/feedback heating capture the key thermal physics 
of cool-core clusters. 



5 ADDITIONAL ASTROPHYSICAL 
IMPLICATIONS 

Our work on local thermal instability in global thermal equi- 
librium is relevant to several other puzzles related to the hot 
gas in massive halos. Here we discuss three applications in 
detail. 

5.1 AGN Fueling: Hot or Cold Gas? 

For massive elliptical galaxies, Allen ct al. (2006) argued that 
the observed jet power (calculated from estimating the energy 
and inflation timescale for X-ray cavities/radio-bubbles) is 
well correlated with the estimated Bondi accretion rate of 



the hot ICM (see Narayan & Fabian 2011 for models of this 
physics). The estimate of the Bondi accretion rate is quite 
uncertain because it requires extrapolating the density and 
temperature from large radii to the radius of influence of the 
supermassive black hole. Nonetheless, the possible correlation 
of the Bondi accretion rate and the jet power suggests that 
the supermassive black hole may be fed by the hot ICM. In 
our calculations, we indeed find that systems spend most of 
their time accreting primarily via the hot phase (Fig. 2). 

For several reasons, however, it is not clear whether the 
connection between the Bondi accretion rate and the true 
AGN accretion rate is that robust. The Bondi accretion rate is 
cx K^'^^^ (where K is the entropy at the Bondi radius). Thus 
lower entropy implies a larger Bondi accretion rate. However, 
lower entropy also implies smaller tcooi/iff (see Eqs. 8-10 for 
the definition of different important timescales) and thus a 
higher accretion rate via cold gas formed by local TI. Cav- 
agnolo et al. (2008) showed that clusters with the central 
entropy < 30 keV cm^ (or equivalently, tTi/tg < 10; see Eq. 
11) show a large Ha luminosity and a large radio luminos- 
ity (a proxy for jet power), suggesting a connection between 
cold gas, AGN feedback, and the state of the hot gas. Indeed, 
our Tl model predicts that the mass accretion rate, and hence 
feedback, increases dramatically whenever tcooi/iff < 10. By 
contrast, the dependence of AGN accretion and feedback on 
the thermal state of the ICM is not expected to be as abrupt 
if accretion proceeds primarily via hot gas. 

For the specific systems in AUcn ct al. (2006), the cen- 
tral gas densities are > 0.1 cm~^. Thus, these systems are 
in fact prone to forming cold filaments at small radii ac- 
cording to our tcooi/iff criterion. Future observations with 
density /temperature measurements close to the black hole's 
sphere of influence are required to ascertain if the Bondi scal- 
ing of the jet power is robust, or if it is by product of the corre- 
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profile of a given halo is isothermalJ In this case 
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Figure 15. The minimum of the ratio for the fiducial run 

and the feedback heating simulations sampled every 50 Myr as a 
function of time (txi/tfi is evaluated from angle averaged radial 
profiles of plasma with T > 0.1 keV; the results do not depend 
sensitively on the cutoff temperature as long as it is < 1/3 the 
ambient temperature). Filaments condense out of the hot phase 
whenever tTl/iff < 10 for both idealized and feedback heating runs. 
This can be seen by comparing this plot with the mass accretion 
rate as a function of time (Figs. 1, 11, & 13). The spike in the mass 
accretion rate corresponds to infalling cold filaments; it is preceded 
by a drop in min(txi/% ) to less than the critical value 10) and 
followed by an increase in min(txi/tff) as the ambient plasma is 
heated and cold matter drops out of the hot phase. 



Lx = 47rrf Ao (T/To 



,1/2 



Jo 



where Ao ~ 10"^'^ erg cm"^ s~^ is the emissivity at To ~ 
10^ K, rs = r2oo/c is the scale radius, c is the concentration 
parameter, and x = r/vs- The integrated X-ray luminosity is 
dominated by the core radius rc, at which the density profile 
transitions from being relatively shallow to steep; the X-ray 



where ric is the gas 



luminosity then scales as Lx oc r^T^^^r 
number density in the core. 

For self-similar evolution of baryons following the dark 
matter (Kaiser 1986) the gas number density is indepen- 
dent of the halo mass; thus, Lx oc r2ooT^^^ oc T^, since 
^200 oc M^/^ and T oc M/r2oo oc M^^^. However, the as- 
sumption of a constant core density irrespective of the halo 
mass is not consistent with our results, especially for low halo 
masses where cooling is very important. In our calculations, 
the X-ray emitting plasma in the core of groups and clusters 
has a lower density for a lower mass halo (e.g., see Fig. 9). 
This is because tcooi/iff 10 irrespective of the halo mass im- 
plies that icooi/iff oc T^^^/nc ~ constant for all halo masses, 
i.e., ric oc T^^^. Assuming for simplicity that the core radius 
follows the self-similar scaling oc M^^^ (which is similar to 
what is observed but not likely to be correct in detail) we find 
that 

X ocn^l ' oc 1 , 

or, more precisely, that the maximum X-ray luminosity at 
a given halo mass is oc . This result is tantalizingly sim- 
ilar to that observed, but a more detailed one-dimensional 
model (e.g., Voit ct al. 2002) is required to quantitatively 
compare our predictions with observations. Such a model will 
be presented in the near future (Sharma et al., in prep.). 



lation between entropy (or, better, tcooi/iff) and the presence 
of cold gas via local TI. 



5.2 Lx - Tx Relation 

The observed relationship between the X-ray luminosity 
and the X-ray emissivity-weighted temperature of groups 
is steeper {Lx oc Tx) than the self-similar scaling predic- 
tion [Lx oc t1; e.g., Evrard & Henry 1991; Bryan 2000 and 
references therein). Here we argue that the non-self similarity 
introduced by the requirement icooi/tft ^ 10 can qualitatively 
account for these observations. 

The total X-ray luminosity produced by a hot halo is 

rr200 

Lx = UeUiAd r, 

Jo 

where r^oo is the virial radius. For simplicity, assume that 
the plasma in groups and clusters emits via thermal free- 
free emission, such that A oc T^^^ , and that the temperature 



5.3 HVCs, TI, & Density of the Galactic Halo 

Much of the gas falling into the Galactic halo is expected to 
virialize because the halo mass ~ 10^^ Mq (e.g., Xue ct al. 
2008) is larger than the critical halo mass required to form a 
virial shock (e.g., Kcrcs ct al. 2009). Local TI, however, can 
generate cold gas if the halo density is high enough. Indeed, 
some fraction of the high velocity clouds (HVCs) of neutral 
Hydrogen observed in the Galactic halo may result from local 
TI (e.g.. Mailer & Bullock 2004). 

Figure 9 shows that for our models of hot halos, a smaller 
halo mass implies that the size of the density core increases 
(relative to the virial radius) and the core density decreases. 
We can generalize this result to Galactic halo masses (with 
temperatures ~ few x 10*' K) in which free-free emission is 
sub-dominant relative to cooling by lines, resulting in an even 
shorter cooling time. Our one-dimensional models of hot ha- 
los based on a tTi/ift threshold (Sharma et al., in prep.) 

^ These assumptions are not strictly correct but do not qualita- 
tively affect our argument. 
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show that the core radius for a Galactic-mass halo is close 
to the virial radius (« 100 kpc) and the core density is 
Tie ~ 10"* cm~"^. This is consistent with the currently poor 
constraints on the properties of the Galactic halo (e.g., An- 
derson & Brcgman 2010). 

Given these parameters, local TI can produce cold 
gas (e.g., HVCs in the Galactic context) at ~ 100 kpc, where 
the ratio tcooi/iff is smallest and close to our threshold for 
forming nonlinear structure. Although the gas condenses at 
large radii it falls inward because of gravity and may be ob- 
served at much smaller radii. Observations in M31, where the 
HVC distances are easier to measure, have found HVCs as far 
out as 50 kpc (and perhaps even farther; Thilker et al. 2004). 
If the HVCs are in pressure balance with the hot halo their 
emissivity is expected to be smaller at large radii where the 
pressures are lower. The distribution of cold g func- 
tion of radius is difficult to predict with certainty because 
of uncertainties in the destruction of the cold filaments by 
conduction, ram pressure stripping, etc. 

Birmcy, Nipoti, & Fratcrnali (2009) argued that the 
Galactic HVCs cannot form due to local TI. As we discuss 
in the Appendix, Birmcy, Nipoti, & Fratcrnali (2009) 's ar- 
guments are appropriate for homogeneous cooling-flows with 
very small amplitude density perturbations. For halos in ap- 
proximate global thermal equilibrium, however, TI can gen- 
erate multi-phase structure (and thus potentially HVCs) if 
iTi/tff <^ 10. Even cooling- flows can lead to multi-phase gas 
if both tTi/ta is small enough and there are large amplitude 
initial density perturbations (see Appendix A for details). We 
discuss the differences between our results and those of Bin- 
ncy, Nipoti, & Fratcrnali (2009) in more detail in Paper I. 

Although cooling-flow models can, under some circum- 
stances, produce multi-phase structure analogous to HVCs, 
we have shown that the presence of approximate global ther- 
mal equilibrium in hot halos makes the formation of multi- 
phase structure far more generic, e.g., even with very small 
density perturbations. This is true even if there are substan- 
tial (order unity) fluctuations in heating and cooling about 
the global equilibrium (see Paper I and section 3.5). 



6 SUMMARY 

The primary source of baryons in galaxies is the smooth accre- 
tion of intergalactic gas into dark matter halos (as opposed to 
mergers; this is particularly true at the current redshift; e.g., 
see Dekel et al. 2009; Keres et al. 2009). For halo masses 
< lO^'^'^ Mq (such that the cooling time is < the free-fall 
time at the virial radius; e.g., Rccs & Ostrikcr 1977) the ac- 
creted baryons remain cold inside the virial radius of the halo. 
For more massive halos, the cooling time at the virial radius 
is longer than the free-fall time and the virial shock creates 
pressure-supported hot plasma that fills the dark matter halo. 
The focus of this paper is on this latter mode of accretion, 
relevant for clusters, groups, and massive galaxies. Although 
the cooling time is longer than the free-fall time at the virial 
radius, the cooling time can be significantly shorter than the 



Hubble-time (and in principle even shorter than the free-fall 
time) at small radii if the gas density profile were to approx- 
imately follow the dark matter halo density profile. Thus, 
cooling has a profound impact on the hot ICM in massive 
halos. 

If the cooling time is shorter than the 'age' of a 
halo (which is approximately the time since the last major 
merger in a realistic cosmological context), a massive cooling- 
fiow will develop absent additional heating. Thus, the obser- 
vational absence of cooling-flows (e.g., Peterson & Fabian 
2006) and the lack of star formation rates of the order of the 
expected cooling rates (e.g., Raffcrty, McNaniara, & Nulsen 
2008) imply that there is an additional heat source roughly 
balancing cooling in the hot halos of groups and clusters. 
Other proposals such as a uniform preheating (e.g., Tozzi 
& Norman 2001) may be important for systems with high 
central entropies and long central cooling times, but cannot 
explain the presence of long-lived cool-core clusters (e.g., Mc- 
Carthy et al. 2008). 

While there are many potential sources of energy that 
can maintain approximate thermal equilibrium, the particu- 
larly attractive feature of heating by star formation and/or 
an AGN is that it can be globally thermally stable (e.g., see 
Guo ct al. 2008 for a specific realization of this): the in- 
creasing feedback due to the cooling and infiow of gas pro- 
vides a thermostat which prevents the ICM core from run- 
away cooling. In this paper we have addressed a previously 
under-appreciated aspect of this scenario: even if heating can 
balance cooling in a globally stable manner, and even if the 
heating can be distributed over the full volume of the ICM (a 
difficult problem; Vcrnalco & Reynolds 2006; Scannapieco & 
Briiggon 2008; Dubois ct al. 2010; Gaspari et al. 2011), 
the cooling of the ICM is very likely to generate local thermal 
instability (TI). 

We have addressed a number of key science questions as- 
sociated with local TI in hot halos. The first is the condition 
under which local TI in a globally stable ICM produces a 
multi-phase medium with cold gas co-spatial with the hot 
ICM. We find that the criterion for TI to produce multi- 
phase gas in global spherical numerical simulations is that the 
thermal instability growth time tTi satisfy tTi < 10 tg (see 
Fig. 8 & the left panel of Fig. 1), where ta is the local free-fall 
time (note that this criterion is derived from global simula- 
tions but it applies at each location within a hot halo; see Eqs. 
8-10 for the definition of these timescales). When tTi > Wtg, 
the TI leads to density and temperature fluctuations but they 
do not become nonlinear and so multi-phase gas cannot gener- 
ically be produced from arbitrary initial conditions. 

In Paper I, we presented a closely related study of local 
TI in globally stable systems using analytic theory and local 
Cartesian simulations. A more comprehensive discussion of 
the key physics is given there. In the local simulations, the 
criterion for TI to produce multi-phase gas is that Iti <, tft 
(rather than txi < 10 tg as we find here). Thus it is easier 
for cold gas to condense out of the ICM in spherical geom- 
etry than in Cartesian geometry. Our interpretation is that 
this difference is because in spherical geometry a slightly over- 
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dense fluid element is compressed as it falls; this increases the 
density and reduces the cooling time, aiding the development 
of the TI.*^ Our criterion for the production of multi-phase 
gas via TI does not depend sensitively on how we implement 
the uncertain heating of the ICM. We have studied two differ- 
ent models: an 'idealized' heating model in which the heating 
in a given radial shell is equal to the average cooling within 
that shell (Eq. 6) and a 'feedback' heating model in which 
the heating is set by the mass accretion rate at the small- 
est radii (Eq. 7). Both heating models produce similar results 
overall and quantitatively similar criterion for the production 
of multi-phase gas via TI. 

Figure 12 in Paper I shows the observed radial profiles of 
icooi/tfi for a number of clusters, some of which show spatially 
extended Ha emission and some of which do not. The tran- 
sition between clusters with and without this observational 
signature of intracluster cold gas occurs at tcooi/tff ~ 5 — 10, 
consistent with our results on the origin of multi-phase gas 
via TI in the ICM. 

The formation of multi-phase gas via TI in hot halos 
makes a number of other predictions about the observed prop- 
erties of the ICM that are in reasonable agreement with obser- 
vations. We highlight several of the key ones here. First, the 
net cooling rate to low temperatures (and the mass accretion 
rate at small radii) for models in global thermal balance is 
10^~^ times smaller than for cooling-fiow models, consistent 
with observations (e.g., Peterson ct al. 2003; O'Dca ct al. 
2008). The net cooling rate is the largest when tTi/tff < 10 
and cold gas is condensing out of the hot phase. Even for mod- 
els with tTi ~ 10 iff , however, which is typical of cool-core 
clusters, the net cooling rate is a factor of ~ 100 smaller than 
the cooling-flow value for that cluster (e.g.. Fig. 3). This is a 
nontrivial result that is a consequence of both global thermal 
balance and the physics of local TI in gravitationally strati- 
fied plasmas. A corollary of this result is that in a given halo 
the observed X-ray emission is produced almost entirely by 
gas within a factor of a few of the virial temperature (Fig. 10). 
This is consistent with observational results from X-ray spec- 
troscopy (e.g., Peterson et al. 2003; Sanders et al. 2010). Pre- 
vious authors have searched for heating mechanisms that can 
self-adjust to reproduce this narrow range in temperatures in 
a given halo (e.g., Kim, El-Zant, & Kamionkowski 2005). Our 
results suggest that these observations are instead a natural 
consequence of the approximate isothermality of halo poten- 
tials, together with the physics of local TI for tTi > 5 — 10 tff , 
which results in only a very small fraction of the gas at any 
time cooling to sub- virial temperatures. 

A second key consequence of the formation of multi-phase 
gas via TI is related to how the global thermal balance of hot 
halos is maintained in the first place. Figure 1 shows that 
the accretion rate at small radii increases dramatically, by 
up to ~ 3 orders of magnitude, when tTi <^ 10 te and cold 

* In addition to spherical compression, the background entropy 
stratification (e.g., Paper I) and the size distribution of cool 
'blobs' (e.g., Pizzolato & Soker 2005) may also afi'ect the formation 
of cold filaments. 



gas condenses out of the hot ICM. In any model in which 
the heating of the hot halo is related to this inflow rate (e.g., 
heating via an AGN or star formation), the heating of the 
ICM will increase dramatically if t^i < 10 iff . This will in turn 
drive tTi/iff 10. As a result hot halos will tend to self-adjust 
such that min(tTi/t£F) ^ 10. Both our 'idealized' heating and 
'feedback' heating models demonstrate this explicitly (Figs. 8 
& 15). We suggest that this physics is responsible for the 
signiflcant fraction of clusters that have central cooling times 
~ 3 X 10* yrs (e.g., Cavagnolo ct al. 2009), which corresponds 
to a minimum value of fcooi ~ 10 tft in cluster cores (see 
Fig. 12 of Paper I). 

The net heating efficiency required to maintain approx- 
imate global thermal equilibrium in clusters is ~ 1 — 3 x 
10"'' (we define the efficiency as the heating integrated over 
the region with a cooling time < 5 Gyr divided by the 
rest mass energy accreted through our inner boundary at 
<; 1 kpc). This efficiency can plausibly be provided by a cen- 
tral AGN. For lower mass halos, the efficiency required for 
thermal equilibrium is lower. Extrapolating to halo masses 
~ 10^*^ M0, an efficiency of < 10"^ would be more than suf- 
ficient for global thermal equilibrium with min(tTi/iff) > 10 
(Sharma et al., in prep.). This is comparable to the efficiency 
of energy injection by core-collapse supernovae, implying that 
supernovae venting from galactic disks may be sufficient to 
maintain global thermal equilibrium in halos like that of the 
Milky Way. 

In this paper we have focused on how cold gas-rich 
filaments condense out of a hot halo; we have not, how- 
ever, addressed important questions about the structure and 
stability of the cold gas. High-resolution numerical simula- 
tions with magnetic fields and anisotropic thermal conduc- 
tion are required to study the survival of the filaments as 
they fall through the hot gas. Arguments based on mixing 
due to Kelvin-Helmholz instability suggest that once a large- 
enough (^ 100 pc) blob cools out of the hot phase it survives 
its passage through the hot medium (e.g., see Murray et al. 
1993; Pizzolato & Soker 2005). Magnetic fields further sta- 
bilize the cold blobs against mixing so we expect the cold 
filaments to survive as they fall toward the center, as is ob- 
served in many cases. 

In section 5 we briefiy summarized the implications of 
our results for other problems related to the astrophysics of 
hot halos, specifically the relative importance of hot and cold 
gas for AGN fueling in clusters (5.1), the non self-similarity of 
the relationship between X-ray luminosity and temperature 
for groups and clusters (5.2), and the possible origin of high 
velocity clouds of neutral Hydrogen in the local group via TI 
in hot halos (5.3). These applications will all be studied in 
more detail in future work. 



ACKNOWLEDGMENTS 

We thank Dusan Keres, Leo Blitz, and Arif Babul for use- 
ful discussions. Support for P. S. was provided by NASA 
through Chandra Postdoctoral Fellowship grant PF8-90054 



© 2011 RAS, MNRAS 000, 1-26 



20 P. Sharma, M. McCourt, E. Quataert, I. J. Parrish 



awarded by the Chandra X-Ray Center, which is operated by 
the Smithsonian Astrophysical Observatory for NASA under 
contract NAS8-03060. M. M., I. P. and E. Q. were supported 
in part by NASA Grant NNX10AC95G, NSF-DOE Grant 
PHY-0812811, and the David and Lucile Packard Founda- 
tion. Computing time was provided by the National Science 
Foundation through the Teragrid resources located at Purdue 
University and NCSA. Some of the runs were carried out on 
Henyey, the theoretical astrophysics computing cluster at UC 
Berkeley supported by NSF AST Grant 0905801. We are also 
thankful for the hospitality of the Kavli Institute for Theo- 
retical Physics (KITP) at UC Santa Barbara, where part of 
this work was performed. 



REFERENCES 

Allen, S. W. et al. 2006, MNRAS, 372, 21 
Anderson, M. E. & Bregman, J. N. 2010, ApJ, 714, 320 
Balbus, S. A. & Soker, N. 1989, ApJ, 341, 611 
Binney, J., Nipoti, C, & Fraternah, F. 2009, MNRAS, 397, 
1804 

Birnboim, Y. & Dekel, A. 2003, MNRAS, 245, 349 

Birzan, L. et al. 2004, ApJ, 607, 800 

Bower, R. G. 2006, MNRAS, 370, 645 

Bryan, G. 2000, ApJL, 544, LI 

Cavagnolo, K. W. et al. 2008, ApJL, L107 

Cavagnolo, K. W. et al., 2009, ApJS, 182, 12 

Ciotti, L. & Ostriker, J. P. 1997, ApJL, 487, L105 

Croton, D. J. et al. 2006, MNRAS, 365, 11 

Dekel, A. & Birnboim, Y. 2008, MNRAS, 383, 119 

Dekel, A. et al. 2009, Nature, 2009, 457, 451 

Dubois, Y., Devriendt, J., Slyz, A., & Teyssier, R. 2010, 

MNRAS, 409, 985 
Evrard, A. E. & Henry, J. P. 1991, ApJ, 383, 95 
Fabian, A. C. 1994, Ann. Rev. Astr. Astr., 32, 277 
Ferland, G. J. et al. 2009, MNRAS, 392, 1475 
Field, G. B. 1965, ApJ, 142, 531 

Gaspari, M., Melioh, C, Brighenti, F., & D'Ercole, A. D. 

2011, MNRAS, 411, 349 
Guo, F., Peng Oh, S., & Ruszkowski, M. 2008, ApJ, 688, 859 
Hayes, J. et al. 2006, ApJ, 165, 188 

Heckman, T. M., Baum, S. A., van Breugel, W. J. M., & 

McCarthy, P. 1989, ApJ, 338, 44 
Hicks, A. K., Mushotzky, R., & Donahue, M. 2010, ApJ, 719, 

1844 

Hoyle, F. 1953, ApJ, 118, 513 

Hu, E. M., Cowie, L. L., & Wang, Z. 1985, ApJS, 59, 447 
Johnstone, R. M., Allen, S. W., Fabian, A. C, & Sanders, 

J. S. 2002, MNRAS, 336, 299 
Kaiser, N. 1986, MNRAS, 1986, 222, 323 
Kaufmann, T. et al. 2009, MNRAS, 396, 191 
Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, 

MNRAS, 363, 2 
Keres, D., Katz, N., Fardal, M., Dave, & Weinberg, D. H., 

& Dave, R. 2009, MNRAS, 395, 160 



Kim, W., El-Zant, A., & Kamionkowski, M. 2005, ApJ, 632, 
157 

Kim, W. & Narayan, R. 2003, ApJ, 596, 889 

Mailer, A. H. & Bullock, J. S. 2004, MNRAS, 355, 694 

McCarthy, L G., Babul, A., Bower, R. G., & Balogh, M. L. 

2008, MNRAS, 386, 1309 
McCourt, M., Sharma, P., Quataert, E., & Parrish, I. J. 

2011, arXiv:1105.2563 (Paper I) 
McDonald, M., Veilleux, S., Rupke, D. S. N., & Mushotzky, 

R. 2010, ApJ, 721, 1262 
McDonald, M., Veilleux, S., Rupke, D. S. N., Mushotzky, R., 

& Reynolds, C. S. 2011, ApJ, 734, 95 
McKee, C. F. & Ostriker, J. P. 1977, ApJ, 218, 148 
Molendi, S. & Pizzolato, F. 2001, ApJ, 560, 194 
Murray, S. D., White, S. D. M., Blondin, J. M., & Lin, D. 

N. C. 1993, ApJ, 407, 588 
Narayan, R., & Fabian, A. C. 2011, arXiv:1105.0594 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 

490, 493 

Nulsen, P. E. J. 1985, MNRAS, 221, 377 
O'Dea, C. P. et al. 2008, ApJ, 681, 1035 
Parrish, L J., Quataert, E., & Sharma, P. 2009, ApJ, 703, 
96 

Peek, J. E. G., Putman, M. E., & Sommer-Larsen, J. 2008, 
674, 227 

Peterson, J. R. et al. 2003, ApJ, 590, 207 
Peterson, J. R. & Fabian, A. C. 2006, Phys. Rep., 427, 1 
Pizzolato, F. & Soker, N. 2005, ApJ, 632, 821 
Ponman, T. J., Cannon, D. B., & Navarro, J. F. 1999, Na- 
ture, 397, 135 
Pratt, G. W. et al. 2010, A&A, 511, A85 
Rafferty, D. A., McNamara, B. R., & Nulsen, P. E. J. 2008, 

ApJ, 687, 899 
Rees, M. J. & Ostriker, J. P. 1977, MNRAS, 179, 541 
Sanders, J. S., Fabian, A. C, Frank, K. A., Peterson, J. R., 

& RusseU, H. R. 2010, MNRAS, 402, 127 
Scannapieco, E. & Briiggen, M. 2008, ApJ, 686, 927 
Sharma, P., Chandran, B. D. G., Quataert, E., & Parrish, L 

J. 2009, ApJ, 699, 348 
Sharma, P., Parrish, L J., & Quataert, E. 2010, ApJ, 720, 
652 

Soker, N. & David, L. P. 2003, ApJ, 589, 770 

Sommer-Larsen, J. 2006, ApJL, 644, LI 

Stone, J. M. & Norman, M. L. 1992, ApJS, 80, 753 

Sun, M. 2009, ApJ, 704, 1586 

Sun, M. et al. 2009, ApJ, 693, 1142 

Sutherland, R. S. & Dopita, M. A. 1993, ApJS, 88, 253 

Tamura, T. et al. 2001, a, 365, L87 

Thilker, D. A. et al. 2004, ApJL, 601, L39 

Tozzi, P. & Norman, C. 2001, ApJ, 546, 63 

Vernaleo, J. C. & Reynolds, C. S. 2006, ApJ, 645, 83 

Voit, G. M., Bryan, G. L., Balogh, M. L., & Bower, R. G. 

2002, ApJ, 576, 601 
Wakker, B. P. & van Woerden, H. 1997, Ann. Rev. Astr. 

Astr., 1997, 35, 217 
Werner, N. et al. 2010, MNRAS, 407, 2063 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 



© 2011 RAS, MNRAS 000, 1-26 



Feedback Regulation of Hot Halos 



Xue, X. X. et al. 2008, ApJ, 684, 1143 

Zakamska, N. L. & Narayan, R. 2003, ApJ, 582, 162 



© 2011 RAS, MNRAS 000, 1-26 



22 P. Sharma, M. McCourt, E. Quataert, I. J. Parrish 



Table 1. Runs with Idealized Heating 



Label* 






a' 


till A 


A>** 


^-^0,01-0.1 




^'-'>n.i 




^^^0.01-0,1 






ett = 


Q + /Mc^ 


MCI 


1.91 


209.7 


1.4 


2.9 (2.3) 


7.3 (70.3) 


0.02 (1) 


12.5 (470.6) 


9.4 X 10" (9 > 


; 10") 


4.3 X 10 (3.8 X 10 ) 


4.3 


X 10 (8.1 X 10 ) 


8.2 


X 10-" 


MClO 


19.1 


209.7 


1.4 


32.5 (12.5) 


2.6 (41.4) 


(0.5) 


(348.3) 


8.8 X 10" (8.9 


X 10") 


(3.5 X 10^) 




(6.7 X 10') 


3.7 


X 10-" 


MC30 


57.2 


209.7 


1.4 


111.4 (32.5) 


0.8 (21.1) 


(1.2) 


(169) 


8 X 10" (8.4 ; 


< 10') 


(1.5 X lO**) 




(4 X 10") 


4.3 


X 10-" 


MClOO 


190.6 


285.9 


0.8 


477 (128.5) 


0.3 (0.4) 


(0) 


0(0) 


6 X 10" (6.3 > 


; 10") 


0(0) 




0(0) 







MC300 


,571.8 


152.5 


0.8 


1495 (191) 


0.2 (0.5) 


(0) 


0(0) 


5.1 X 10" (5.2 


X 10") 


0(0) 




0(0) 







CI 


1 


110 


1.4 


1.6 (1.4) 


2.8 (13.7) 


0.03 (3.7) 


9.1 (286.6) 


3.5 X 10" (3.4 


X 10") 


4.3 X 10'' (4.3 X 10") 


3.5 


X 10" (5.6 X 10') 


4.5 


X 10-* 


CIO* 


10 


110 


1.4 


20.2 (8.2) 


1.2 (11.4) 


0.003 (3.04) 


2.3 (230.3) 


3.3 xlO" (3.3 


X 10") 


2.6 X 10' (3.6 X 10") 


1.3 


X 10" (4.9 X 10') 


10"^ 


C30 


30 


110 


1.4 


71 (21.6) 


0.4 (7.7) 


(1.2) 


(144.6) 


3x10" (3.2 X 


10") 


(2.5 X 10") 




(4 X 10') 


5 > 


: 10-" 


ClOO 


100 


150 


0.8 


314.4 (85) 


0.1 (0.3) 


(0) 


0(0) 


2.3x10" (2.5 ) 


< 10") 


0(0) 




0(0) 


5.4 


X 10-= 


C300 


300 


80 


0.8 


1001 (126.4) 


0.1 (0.2) 


(0) 


0(0) 


2 X 10" (2 X 


10") 


0(0) 




0(0) 







MGl 


0.5 


57.7 


1.4 


0.52 (0.46) 


0.4 (2.7) 


0.02 (0.3) 


16.2 (130.4) 


8.9 xlO" (8.6 


X 10") 


3.9 X lO'' (3.2 X 10^) 


1.3 


X 10' (3.3 X 10') 


7 > 


: 10-5 


MGIO 


5.3 


57.7 


1.4 


9.9 (3.8) 


0.3 (2.4) 


0.02 (0.2) 


4.7 (113.9) 


8.6 X 10" (8.4 X 10") 


2 X 10' (2.8 X 10') 


2.91 


X 10" (2.9 X 10') 


1.8 


X 10-* 


MG30 


15.7 


57.7 


1.4 


37.1 (10.7) 


0.2 (1.9) 


0(0.1) 


(89.5) 


7.9 X 10" (8.2 


X 10") 


(2.2 X 10') 
(1.7 X 10') 




(2.4 X 10') 


3.2 


X 10-" 


MGIOO 


52.5 


78.7 


0.8 


169.1 (45.8) 


0.04 (0.6) 


(0.02) 


(27.4) 


6 X 10" (7.3 > 


: 10") 




(1.4 X 10') 


1.2 


X 10-" 


MG300 


157.4 


42 


0.8 


557.4 (68.7) 


0.02 (0.02) 


(0) 


0(0) 


5.2 X 10" (5.9 


X 10") 


0(0) 




0(0) 







G5 


1.1 


23.7 


1.4 


1.8 (0.9) 


0.07 (0.5) 


0.004 (0.07) 


8.6 (73.7) 


3.2 X 10" (3.1 


X 10") 


2 X 10' (3.8 X 10') 


8.8 


X 10" (2.5 X 10') 


4 > 


: 10--'^ 


G50 


10.8 


23.7 


1.4 


39.7 (10.4) 


0.05 (0.4) 


0.003 (0.05) 


2.3 (48.9) 


2.7 X 10" (3 > 


; 10") 


1.2 X 10' (2.3 X 10') 


2 > 


< 10" (1.6 X 10') 


7.4 


X 10-" 


G75 


16.2 


23.7 


1.4 


65.8 (16.2) 


0.03 (0.13) 


(0.4) 


(40.9) 


2.6 X 10" (2.9 


X 10") 


(1.1 X 10') 




(1.3 X 10') 


4.7 


X 10-" 


G150 


32.3 


17.2 


0.8 


158.6 (26.4) 


0.01 (0.12) 


(0.2) 


(41.6) 


2.3 X 10" (2.9 


X 10") 


(9.7 X 10") 




(1.5 X 10') 


5.3 


X 10-" 


G300 


64.6 


17.2 


0.8 


350.3 (42.3) 


0.006 (0.05) 


(0.04) 


(3.9) 


2 X 10" (2.3 > 


: 10") 


(6 X 10') 




(7.8 X 10') 


1.5 


X 10-* 



* CIO is the fiducial run. 

* The entropy K (Tj^cv/'^e''^) is initialized via Eq. .5. 

t CI is the label for a cluster (Mq = 3.8 X 10^"* Mq) with Kq = 1 keV cm^. MCI stands for a massive cluster with Mq = 10^^ M0 and 
a core entropy Kq scaled according to self-similarity, K <x M^^^ , from the cluster model with Kq = 1 keV cm^. MG stands for a massive 
group with Afo = 10^* Mq and G stands for a group with Mq = 3.8 X 10^^ Mq . Initial conditions for all halo masses are scaled with Mq 
according to self-similarity. Self-similarity is broken with time because of cooling and feedback. 

X initial value of the ratio iTi/iff evaluated at the inner boundary (value in the parenthesis is the minimum value over all radii), tff at the 
inner boundary is ~ 0.014 Gyr for all cases. 

ft Af>o.i represents the average (over the whole run for 5 Gyr) mass accretion rate (in Mq yr~^) through the inner boundary for plasma 
with T > 0.1 keV, A/0.01— O.i for 0-01 keV < T < 0.1 keV, and A/<o.oi for T < 0.01 keV. Values in parentheses correspond to 'cooling-flow' 
models with no external heating. Similarly, A/>o.i is the gas mass (in Mq) in the hot phase, and likewise for A/o. 01-0.1 and Af<o.oi; values 
in parentheses correspond to cooling-flow models. 

tl e represents the time-averaged efficiency required to balance cooling; e = j Mc? , where Af is the time-averaged total mass accretion 
rate through the inner boundary, and = J neniA{T)dV and the integral is over all radii where the angle-averaged cooling time is 
< 5 Gyr. 
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Table 2. Additional Cluster Runs 



LabelT 


Nr X Ng X 


M>o.i 


A^O.Ol-O.l 


M<o.oi 


M>o.i 


Mo. 01-0.1 


M<o.oi 


e = Q+/Mc^ 


ClO-pl 


512 X 256 X 1 


1.1 


0.01 


1.7 


3.3 X 10^^ 


1.7 x 10^ 


6.2 X 10^ 


1.3 X lO"'' 


ClO-plO 


512 X 256 X 1 


1.7 


0.001 


0.24 


3.3 X IQi^ 


3.5 X 10'' 


4.8 X 10^ 


2 X 10"'^ 


ClO-nzl 


512 X 256 X 1 


0.3 


0009 


0.089 


4.2 X 10^-"^ 


6.1 X 10^ 


1.1 X 10^ 


3.8 X 10~^ 


C10-nz4 


512 X 256 X 1 


1.1 


0.004 


2.6 


2.7 X 10^^ 


2.9 X 10^ 


1.9 X 10^ 


4.1 X 10~^ 


C10-nz8 


512 X 256 X 1 


1.3 


0.02 


2.5 


3.1 X 10^2 


2.6 X 10^ 


108 


7.1 X 10"^ 


C10-nzl6 


512 X 256 X 1 


1.4 


0.001 


0.44 


3.3 X 1012 


1.3 X 10^5 


1.2 X lO'^ 


2 X 10-3 


C10-3D 


256 X 128 X 32 


0.77 


0.018 


0.95 


3.5 X 10^2 


9.1 X 10^ 


1.6 X 10^ 


2.4 X 10-3 


ClO-hlf 


256 X 128 X 1 


1.3 


0.004 


1.1 


3.5 X 10^2 


4.8 X 10^5 


3.5 X lO'^ 


3 X 10-3 


CIO 


512 X 256 X 1 


1.2 


0.003 


2.3 


3.3 xl0i2 


2.6 X 10^ 


1.3 X 10^ 


10-3 


ClO-dbl 


1024 X 512 X 1 


1.4 


0.002 


0.85 


2.6 X 10^2 


5.4 X 10^5 


2.6 X 10^ 


7 X 10-^ 


C30-pl 


512 X 256 X 1 


0.38 








3 X 10^2 








4.9 X 10-3 


C30-pl0 


512 X 256 X 1 


0.58 








2.9 X 10^2 








2.8 X 10-3 


C30-3D 


256 X 128 X 32 


0.24 








3.2 xl0i2 








8.6 X 10-3 


C30-hlf 


256 X 128 X 1 


0.28 








3.2 xl0i2 








7.6 X 10-3 


C30 


512 X 256 X 1 


0.4 








3x10^2 








5 X 10-3 


C30-dbl 


1024 X 512 X 1 


0.36 








3 x 10^2 








1.5 X 10-3 



fClO-plO is a cluster run with random perturbations ~ 10 (drawn from a uniform distribution) in the heating rate on top of the balance 
of cooling and heating in each shell; the random perturbations are changed every 100 Myr. C10-nz4 corresponds to a run where the 
balance of heating and cooling is applied over 4 logarithmically spaced radial zones (i.e., H/r = 1/4, where H is the scale over which 
heating is applied uniformly). C10-3D is a three-dimensional run with the same parameters as the fiducial run CIO. ClO-dbl (ClO-hlf) has 
double (half) the fiducial resolution. All remaining quantities are defined in the same way as Table 1. 



Table 3. Cluster and Group Runs with Feedback Heating 



Labelt 






M>0.1 


A^O. 01-0.1 


M<o.oi 


M>o.i 


-Mo. 01-0.1 


M<o.oi 


e = Q+/Mc2 


C10-fb-7em3 


7 X 10" 


■3 


0.15 


0.0017 


0.6 


2.84 X 10^2 


8.28 X 10** 


1.07 X 10'^ 


1.4 X 10" 


-3 


C10-fb-4em3 


4 X 10" 


■3 


0.46 


0.006 


2.74 


2.74 X 10^2 


2.02 X 10^ 


1.3 X 10^ 


5.8 X 10" 


-4 


C10-fb-2em3 


2 X 10" 


■3 


0.77 


0.0055 


2.34 


3.18 X 10^2 


1.17 X 10^ 


1.05 X 10^ 


7.2 X 10" 


-4 


C10-fb-7em4 


7 X 10" 


-4 


1.67 


0.019 


4.9 


3.22 X 10^2 


3.5 X 10^ 


2 X 10^ 


3.7 X 10" 


-4 


G50-fb-4em3 


4 X 10" 


■3 


0.022 


0.00042 


0.079 


1.8 X 10^1 


1.88 X 10^ 


5.79 X 10*5 


2.9 X 10" 


-4 


G50-fb-2em3 


2 X 10" 


■3 


0.029 


0.00041 


0.19 


1.76 X 10" 


1.25 X 10^ 


7.24 X 10^ 


1.1 X 10" 


-4 


G50-fb-7em4 


7 X 10" 


-4 


0.026 


0.00046 


0.22 


2.54 X 10" 


4.3 X 10^ 


1.09 X 10'^ 


5.3 X 10" 


-4 


G50-fb-7em5 


7 X 10" 


■5 


0.051 


0.0058 


3.23 


2.99 X 10" 


7.83 X 10^ 


2.18 X 10^ 


7.5 X 10" 


-5 



tC10-fb-7em3 stands for a run with the same initial profile as CIO (the fiducial run) but using a feedback heating model with ef^ = 
7 X 10-=^ (see Eq. 7). The group runs (with the same initial conditions as G50) with feedback heating are also included. All remaining 
quantities are defined in the same way as Table 1. 
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APPENDIX A: INHOMOGENEOUS COOLING 
FLOWS 

In this Appendix we discuss the nature of hot halos in the 
absence of external heating — in this case the halo plasma set- 
tles into a cooling-flow. We consider both steady homogeneous 
cooling-flows and cooling-flows with large initial density per- 
turbations that can develop multi-phase structure. Although 
observations imply that hot, virialized halos are in approx- 
imate thermal equilibrium, the results in this Appendix are 
useful for comparison to other results in the literature and for 
comparison to the simulations with heating described in the 
main text. 



Al Steady Homogeneous Cooling Flows 

In a steady (d/dt = 0) cooling-flow the equations of conser- 
vation of mass, momentum, and energy can be expressed as 



47rr pv = M, 

1 d\nv 1 d\np 2 
"l^dh^" tCrfhi^'^tf 
1 dln(p/pT) 



tu 



din r 



1 

f-cool 



(Al) 
(A2) 

(A3) 



where i; (> 0) is the radial infall velocity, M is the steady mass 
inflow rate, tin = r/v is the infall time, tsnd = 'r/{p/pY^^ 
is the isothermal sound-crossing time, ts is the free-fall 
time (Eq. 10), and tcooi is the cooling time (Eq. 9). We use 
Eqs. A1-A3 to interpret our cooling- flow results discussed 
shortly. 

For our homogeneous cooling-flow simulations we start 
with an initial hydrostatic profile identical to the fiducial sim- 
ulation described in the main text (CIO), but without den- 
sity/temperature perturbations. As in our simulations with 
heating, we use cosmologically realistic infiow boundary con- 
ditions. Since the cooling time is the shortest at the cen- 
ter (but tcooi S> isnd ~ iff in the initial state), gas starts 
to cool and fiow in subsonically (eq A3). The left panel of 
Figure Al shows the key timescales in the resulting steady 
cooling-flow. As the gas cools and flows in to small radii, the 
cooling time becomes less than the free-fall time at r < 5 kpc. 
The gas then abruptly cools to low temperatures and the in- 
flow becomes supersonic; this corresponds to the rapid in- 
crease (decrease) in the sound-crossing (cooling) time in the 
left panel of Figure Al at ~ 1.5 kpc. For a given gravitational 
potential, the radius at which the steady cooling-flow becomes 
supersonic is determined primarily by the density at the outer 
boundary; for a larger density at large radii the sonic radius 
moves out, while for a lower density at large radii, the sonic 
radius decreases and the flow may never become supersonic. 



A2 Multiphase Gas in a Cooling Flow 

Small amplitude density perturbations do not grow signifi- 
cantly via thermal instability in a cooling-fiow because the 
perturbations are advected inwards on the same timescale 



that the thermal instability would nominally develop (Bal- 
bus & Sokcr 1989). This is a critical difference between 
cooling-flows and the plasmas in approximate global thermal 
equilibrium considered in the main text. One implication of 
this result is that cooling-flows with small amplitude density 
perturbations cannot, in general, develop multi-phase struc- 
ture. With sufficiently large density perturbations, however, 
a cooling-fiow can produce multi-phase gas (e.g., Pizzolato & 
Sokcr 2005). The resulting structure is very different from 
the homogeneous cooling-flows described in the previous sec- 
tion. A rough criterion for whether a cooling-fiow becomes 
multi-phase can be derived as follows. Multiphase gas will 
result if the effective cooling rate of an overdense blob (the 
cooling rate of the blob minus the cooling rate of the back- 
ground cooling-fiow) is larger than the effective free-fall rate. 
The effective free-fall rate is 

1/2 

(A4) 



iff,cff 



1 



where 5 = \5p/ p\ ^ \5T/T\ is the overdensity of the blob 
relative to the ambient density. The effective cooling rate of 
an overdense blob is 

1 1 1 

(A5) 



^cooljCff icool,b icool,a 

where icooi.b is the cooling rate of the blob and tcooi.a is the 
cooling rate of the ambient medium. For small overdensities, 
and in the isobaric limit, the effective cooling rate can be 
simplified to 



d 



ainT 



-1 

^cool 



5_5_ 

3fTl' 



where t^i is the isobaric TI timescale given by Eq. 8. The 
effective free-fall rate for small overdensities is l/tfi,cff ^ 
5^^'^ /ts, so that a cooling fiow can become multi-phase if 



1/2 



iff 



(A6) 



(again assuming 5 <^ 1). Note also that a cooling-flow can 
become multi-phase only if the cooling curve is thermally un- 
stable in the isobaric regime; i.e., if dluA/dlnT < 2; when 
this condition is not met, the slightly overdense blobs are 
heated relative to the background and hence cannot become 
denser/cooler. 

The right panel of Figure Al shows the angle averaged 
timescales in the quasi-steady state for a cooling-fiow sim- 
ulation with density perturbations identical to the fiducial 
run ((5max ~ 0.3); these timescales are calculated using the 
properties of the hot plasma only (T > 0.1 keV). The sound 
crossing time for the hot plasma is < the free-fall time at all 
radii, implying that the hot plasma is in approximate hydro- 
static equilibrium at all radii. Within r <7 kpc, however, the 
cooling time in the hot plasma, which is of order the infall 
time (see Eq. A3), is also of order the free-fall time. While 
the conservation laws of momentum and energy (eqs. A2 & 
A3) are roughly valid when applied solely to the hot phase, 
the mass conservation equation (Eq. Al) is not. Specifically, 
mass is not conserved for the hot plasma within r <^ 7 kpc 
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Figure Al. Key timescales in a steady liomogeneous cooling- flow {left panel) and a quasi-steady multi-pliase cooling- flow with large density 
perturbations [right panel; see eqs. A1-A3 for tlie precise definition of tlie different timescales). Only the plasma in the hot (T > 0.1 keV) 
phase has been used to calculate the relevant timescales for the multi-phase cooling-flow simulation. In the right panel the sound crossing 
time (tsnd) and the infall time (tin) are only shown for our fiducial density perturbations (5max = 0.3). 



because mass "drops out" via thermal instability to the cold 
phase: the mass accretion rate in the hot plasma thus de- 
creases with a decreasing radius for r <7 kpc. This is shown 
explicitly in Figure A2 which shows the average mass accre- 
tion rate in the hot phase (T > 0.1 keV) as a function of 
radius for the cooling-flow simulations with different initial 
density perturbations. 

The right panel of Figure Al also shows the cooling 
timescale in the hot plasma in the quasi-steady state for the 
cooling flow simulations initialized with density perturbations 
3 times bigger (5max = 1) and 10 times smaller ((5max = 0.03) 
than our flducial case. The ratio of the cooling time to the 
free-fall time at the radius where the cooling-flow becomes 
multi-phase increases with larger density perturbations and 
the multi-phase cooling flow sets in at larger radii (see Eq. 
A6). 

Figure A3 shows the electron number density in the in- 
ner 10 kpc for the quasi-steady state of cooling-flow simula- 
tions with different initial density perturbations. Figure A3 
shows the qualitative difference between a cooling-flow with 
and without significant initial density perturbations: cooling- 
flows without significant density perturbations do not show 
dense multi-phase gas (consistent with the analytic arguments 
of Balbus & Soker 1989); instead, the plasma just monolithi- 
cally cools to low temperatures at radii where the cooling time 
is less than the free-fall time (see Fig. Al). Figure A3 also 
visually demonstrates that the multi-phase structure that de- 
velops in cooling-flows with large initial density perturbations 
depends strongly on the magnitude of the density perturba- 
tions: the models with larger density perturbations develop 
more pronounced multi-phase structure out to larger radii. 
This is a critical difference between cooling-flow models and 
models in which the plasma in hot halos is in approximate 
global thermal equilibrium. In the latter, the thermal insta- 
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Figure A2. Mass accretion rate in the hot phase (T > 0.1 keV) as 
a function of radius averaged from 4 to 5 Gyr. The total accretion 
rate is relatively independent of both radius and the magnitude 
of the density perturbations (parameterized by (5max; see Fig Al). 
However, the mass accretion rate in the hot phase decreases at 
small radii by an amount that depends on the amplitude of the 
density perturbations; larger density perturbations imply more gas 
cooling to low temperatures and thus a lower mass accretion rate 
in the hot phase. 



bility grows exponentially and the multi-phase structure that 
develops is essentially independent of the initial density per- 
turbations. Thus the development of multi-phase structure 
is far more generic and robust if the virialized plasma is in 
approximate global thermal equilibrium. 

Cooling-flow simulations are quite sensitive to the bound- 
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Figure A3. Snapshots of the electron number density in the inner 
10 kpc for the coohng-fiow simulations with different fractional 
density perturbations S = \5p/p\. The simulations are initialized 
with a distribution of 5 identical to our fiducial run but rescaled 
such that <5 < <5niax- This Figure demonstrates that a cooling-flow 
becomes multi-phase only when the cooling time is short and the 
density perturbations are sufficiently large (see Eq. A6). The color 
scale for the electron density goes from 10"^ '^^ to 10 cm~^ but 
the cold (~ 10** K) gas in <5 ^ simulations is actually much 
denser {~ 100 - 1000 cm-^). 

ary condition at the outer radius because the free-fall time at 
the outer boundary is shorter than the Hubble time. A cos- 
mologically realistic boundary condition should allow mass 
inflow at the outer boundary. Cooling-flow simulations which 
do not allow inflow, e.g., the isolated halo simulations of Kauf- 
niann ct al. (2009), evolve qualitatively differently from sim- 
ulations with inflow such as ours. Without mass supply at 
the outer boundary all the gas with cooling time shorter than 
the age cools and accumulates at the center. Multiphase gas 
in such simulations occurs only for a short time comparable 
to the cooling time. In contrast, our cooling-flow simulations 
which allow inflow at the outer boundary reach a quasi-steady 
state with multiphase gas if density perturbations are large 
enough. 



